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Part  I 


Chapter  I:  Introduction 

The  Relativistic  Vlasov-Maxwell  (RVM)  system  is  a  nonlinear  system 
of  first  order  partial  differential  equations  that  models  the  time 
evolution  of  a  collisionless  plasma,  e.g.  a  high  temperature,  low  density 
ionized  gas.  Numerical  computation  of  solutions  of  this  system  is 
prohibitively  expensive  in  part  because  of  the  six-dimensional  phase 
space  for  the  Vlasov  density  function.  For  computational  feasibility,  we 
consider  a  version  (RVM)  in  which  the  Vlasov  density  f  depends  on  one 
spatial  variable,  x,  and  two  momentum  variables,  v-i  and 
V2.  This  is  the  simplest  version  of  the  problem  which  retains  the 
hyperbolic  structure  of  Maxwell's  Equations  and  for  which  there  is  a 
nontrivial  magnetic  field.  We  treat  the  case  of  a  single  species  of 
particles  with  distribution  function  f,  in  the  presence  or  a  neutralizing 
background  with  density  n(x).  The  electric  field  is  given  by 
E(t,x)=(Ei(t,x),  E2(t,x))  while  the  (scalar)  magnetic  field  is  denoted 
B(t,x).  The  speed  of  light  is  c,  and  we  assume  that  the  rest  mass  and 
charge  of  the  particles  are  both  1.  The  particles  move  under  the  action  of 
their  self-induced  Lorentz  force,  F=(E+c'1vxB),  and  the 
objective  is  to  track  the  simultaneous  evolution  of  the  density  and  the 
fields,  which  satisfy  the  following  Cauchy  Problem: 


(i.i) 


(  RVMJ 


'  A  A 

3  tf  +  Vi3xf  +  (E+C‘1BMv)-Vvf  =  0 

8tEi=-4jtji  dxEi=4rcp 

.3tB=-c3xE2  5tE2=-c3xB-4rcj2 


1 


f 


where  v  =  v(1+<r2|v|2)-1/2  is  the  relativistic  velocity  and 
with  Cauchy  data 

f(x,v,0)  =  fO(x.v)  >  0 
E2(x,0)  -  E°(x) 

B(x,0)  -  B°(x) 

The  charge  and  current  densities,  p  and  j,  are  defined  by 
p(t,x)  =  Jf(t,x,v)dv  -  n(x) 

j(t,x)  =  jvf(t,x,v)dv 

All  data  functions  are  taken  to  be  smooth  and  compactly  supported,  with 
f°  ^  0.  The  background  density,  n(x),  in  addition  to  being  smooth  and  of 
compact  support,  is  neutralizing  in  the  sense  that 

J  p(0,x)dx  =  0 

x 

We  take  Ei(0,x)  =  4 n  Jp(0,y)dy  as  initial  data  for  Ei. 

-oo 

In  [2],  Glassey  and  Schaeffer  proved  global  existence  of  smooth  solutions 
of  this  problem,  their  main  result  being  that  f°  e  Cq(R3)  and 

e£,  B°e  Co(R)  imply  the  existence  of  a  global,  C1  solution  (f,E,B)  that 

satisfies  the  initial  values  f(0,x,v)  =  f°(x,v),  E2(0,x)  =  E^x),  B(0,x)  = 
B°(x)  and  whose  components  f,E,  and  B  are  compactly  supported  V  t  >  0. 
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In  [3],  they  proved  convergence  of  a  particle-in-cell  method  for 
approximating  solutions  of  this  system.  Their  scheme  (as  well  as  other 
typical  schemes  such  as  that  used  in  [1]  ),  is  limited  by  the  constraint 
imposed  by  the  Courant-Friedrich-Lewy  (CFL)  condition  on  the  size  of 
allowable  time  steps,  resulting  in  enormously  expensive  computations. 
CFL  requires  At  <  Ax/c,  since  the  speed  of  light,  c,  is  the  propagation 
speed  for  the  hyperbolic  Maxwell  system  satisfied  by  E2  and  B.  However, 
there  are  regimes  in  which  this  restriction  is  actually  far  too  severe. 
When  the  electromagnetic  quantities  vary  slowly,  it  may  be  possible  to 
replace  Maxwell's  equations  with  simpler  models  that  are  elliptic  in 
nature,  hence  may  lead  to  cheaper  computation  by  allowing  larger  time 
steps.  (  See  for  example  [5].  ) 

In  this  paper,  we  propose  a  modification  of  (RVM),  denoted  (RVM~), 
that  employs  changes  of  this  nature.  The  modifications  are  suggested 
geometrically  by  the  form  of  the  integral  representations  of  the  Maxwell 
fields  E2  and  B  which  result  from  solving  the  1-D  Maxwell  system 

f  3tE2  +  c9xB  +  4n\2  =  0 
|  3tB  +  c9xE  2  =0 

The  Riemann  invariants  are  (E2  +  B)  and  (E2  -  B),  and  we  find 
(1.2a)  E2(t,x)  -  ^  [E2(x-ct)+E2(x+ct)  +  B°(x-ct)-B°(x+ct)]  - 

2*e'[i2]<t,x) 

(1.2b)  B(t.x)  -  2  [E2(x-ct)-E2(x+ct)+B°(x-ct)+B°(x+ct)  ]  - 

2x£C[j2](t,x) 
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t 

where  £±  [j2](t,x)  :  =  J  [  j2(t,x-c(t-x))  ±  j2(x,x+c(t-t))]  dt 

Treating  the  speed  of  light  c  as  a  parameter  and  allowing  it  to  grow  has 
the  effect  of  "flattening  out"  the  cone  integrals  appearing  in  (1.2). 

Formally,  as  c  ->  the  cone  integrals  become  purely  spatial  integrals  and 
we  are  led  to  alternate  field  operators  resulting  in  approximate  fields  E2 
and  B  given  by 

(1.3a)  E2(t,x)  =  ^  [E2(x-ct)+E2(x+ct)+B°(x-ct)-B°(x+ct)]  - 

2*e®[f2](t.x) 

(1.3b)  B(t,x)  -  j  [E2(x-ot)-E2(x+ct)+B°(x-ct)+B0(x+ct)]  - 

2ite?N2](t,x) 

X  00 

where  £±[j  2](t,x)  “  j  Mtyjdy  ±  ^  J  f s(t,y)dy 

-oo  x 

and  j  (t,x)  =  |  vf(t,x,v)dv  is  the  current  density  obtained  from 

R2 

the  density  function  f  of  the  modified  problem.  Note  that  the  same 
Cauchy  data  appears  in  both  (1.2)  and  (1.3),  which  is  natural  given  that  the 
approximation  involved  only  the  integral  operators.  If  we  omit  the  data 

♦v  *SJ 

terms  in  the  expressions  for  E?  and  B,  what  remains,  e.g.  the  integral 
terms,  are  solutions  of  Maxwell’s  equations  as  modified  with  the  Darwin 
and  quasi-electrostatic  approximations.  In  the  one-space,  two-momenta 
case,  this  means  neglecting  the  8tB  and  3tE2  terms  in  (1.1).  However,  if  we 
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make  these  simplifications  of  Maxwell  first,  the  representations  of  the 
solutions  of  the  resulting  system  will  not  contain  these  data  terms.  As 
will  be  shown  later,  the  presence  of  these  terms  results  in  better  (e.g. 
1/c2  vs  1/c)  convergence  of  the  solutions  of  the  two  problems  for  small 
time.  We  propose  the  following  modified  version  of  (RVM): 


(1.4) 


(RVM~) 


9tf  +  Vi3xf  +  ( E+c^BMvJ-V  vf 

0  <  f°(x,v)  6  Cq(R3) 

e£(x)  &  B°(x)  €  Cq(R) 


-  0 


The  data,  f°,  E°  ,  and  B°  are  considered  known.  We  also  assume  a 

smooth,  compactly  supported,  neutralizing  background  density  n(x),  as  in 
(RVM).  Charge  and  current  densities  are  given  by 

p(t,x)  =  J  f  (t,x,v)dv  -  n(x) 

R2 


F  (tx)  =  J  vf  (t,x,v)dv 

R2 


The  electric  and  (scalar)  magnetic  fields  are  given  by 


x 

Ei(t,x)  =  4k  J  p~ (t,y)dy 

-oo 
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E2(t,x)  = 


1 

2 


E^(x-ct)  ••  Eg(x+ct)  +  B°(x-ct)  -  B°(x+ct) 


2k  £®[f  2l(t.x) 


B(t,x)  =  g 


E^(x-ct)  -  Eg(x+ct)  +  B°(x-ct)  +  B°(x+ct) 


-2k  e!tr2](t.x) 

M  M  ^ 

We  assume  existence  of  a  global,  C1  solution  ( f ,  E,  B  ).  The  proof  is 
straightforward  and  similar  to  the  existence  proof  for  the  1-D  Vlasov 
Poisson  system,  once  the  a-priori  bounds  obtained  in  Chapter  II  are  in 
hand. 

In  part  I  of  this  paper,  we  will  show  that  solutions  of  the  modified 
problem  converge  in  a  pointwise  sense  to  solutions  of  the  unmodified 
problem  at  the  asymptotic  rate  of  Me2.  The  main  result  of  part  I  is  the 
following: 

Theorem  :  Let  (f,E,B)  and  (  f,  E,B)  be  the  global,  C1  solutions  of  (RVM) 
and  (RVM~)  respectively,  satisfying  the  same  Cauchy  data 

f°  e  Cq(R3)  and  E°  ,  B°<=  C<J(R). 

For  c  sufficiently  large  compared  to  the  initial  data,  there 
exists  a  non-decreasing  function  D  :[  0,  «>)-»[  o,  °°).  depending 
on  the  initial  data  but  not  on  c,  such  that 
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II  f(0  -  f  (t)  IL  £  ^  for  all  t  >  0 

Part  II  of  the  paper  describes  numerical  experiments  designed  to 
corroborate  this  convergence  rate  and  demonstrate  explicitly  the  value  of 
the  modified  problem  as  an  approximation  to  (RVM). 

As  described  in  Chapter  VI,  a  particle-in-cell  scheme  was  coded  for 
(RVM~)  and  output  was  compared  with  that  obtained  using  the  program 
from  [4].  Besides  demonstrating  the  1/c2  convergence  rate  of  the 
solutions  of  the  two  problems,  the  results  clearly  show  that  the  modified 
problem  does  not  suffer  from  the  CFL  limitation  -  we  achieved  the  same 
accuracy  as  the  Glassey  Schaeffer  scheme  with  significantly  larger 
timesteps. 

Additionally,  to  demonstrate  the  value  of  the  inclusion  of  the  data 
terms  in  the  solution  of  the  modified  problem,  we  compared  output  with  a 
version  of  the  modified  problem  in  which  these  terms  were  omitted,  i.e. 
where  the  expressions  for  E2  and  B  are  solutions  of  the  system  obtained 
by  making  the  Darwin  and  quasi-electrostatic  modifications  of  Maxwell's 
equations.  We  found  that  including  the  data  terms  improved  the  accuracy 
of  the  solutions  with  very  little  increase  in  computation  time.  These 
results  are  discussed  in  more  detail  in  Chapter  VIII. 

Throughout  this  paper,  we  will  use  the  following  conventions  and 
notation: 

•c  is  the  speed  of  light.  We  always  assume  c  >  1. 

•D  will  denote  a  generic  positive  constant  which  depends  on 

the  initial  data,  but  not  on  c. 
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•D(t)  will  denote  a  positive,  non-decreasing  function  with 
domain  [0,°o). 

•Partial  derivatives  will  be  denoted  by  subscripts.  (  for 

a 

example  3X  «  —  ) 

•We  will  frequently  abbbreviate  the  arguments  in  the 

expressions  for  the  characteristics,  writing  for  example 

X(s)  and  V(s)  for  X(s,t,x,v)  and  V(s,t,x,v). 

•We  will  write  f(t)  for  f(t,* , • )  and  f(t,x)  for  f(t,x, * ) . 

•Finally,  we  will  use  the  following  norms: 

For  a  scalar-valued  function  g=g(t,x,v), 

l|g(t)IL  =  sup  {  |g(t,x,v)|  such  that  xe  R  and  vs  R2} 

||  g  IL  =  sup  { |g(t,x,v)|  such  that  t  >  0,  xe  R  and  vs  R2} 

l|gWllc,  -  ItowiL  +  II(»x9>WIL  +  ll(Vvg)(t)IL 

-  llg(t)IL  +  ll(3xg)(t)IL+  ||0V)o)(t)IL+  IIO„g(t)IL 

For  a  vector-valued  function  H  =  H(t,x)  =  (H1(t,x),H2(t,x)), 
l|H(t)IL  =  ||H,(t)IL+  IlHaWlL 
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Chapter  II  :  A-priori  Bounds  for  the  The  Unmodified  Problem 

In  their  proof  of  global  existence  of  smooth  solutions  for  (RVM)  with 
c  =  1,  Glassey  and  Schaeffer  used  a  standard  iteration  scheme, 
convergence  of  the  iterates  being  assured  once  a-priori  C1  bounds  on  f,  E, 
and  B  are  known.  For  the  purposes  of  this  paper,  we  need  corresponding 
bounds  for  the  case  where  c  appears  as  a  parameter  (c  >  1).  In  particular, 
we  must  determine  which  bounds  can  be  taken  to  be  independent  of  c.  We 
assume,  for  any  fixed  c  >  1,  existence  of  a  global-in-time  C1  solution  (f, 
E,  B),  compactly  supported  V  t  >  0.  Proof  of  existence  is  omitted,  since 
the  argument  would  be  essentially  identical  to  that  in  [2],  using  the  a- 
priori  bounds  we  will  find  in  this  chapter,  which  are  obtained  following 
the  methods  in  [2]. 

A.  Density  Estimates 
Define  X(s,t,x,v)  and  V(s,t,x,v)  by 


(2.1a) 


_d_X(s,t,x,v)  =  V^s.t.x.v) 
<  ds 

.X(t,t,x,v)  =  x 


(2.1b) 


d_V(s,t,x,v)  =  E(s,X(s))+c'1B(s,X(s))MV(s) 
<  ds 

.V  (t.t.x.v)  =  v 


Here,  X(s)  and  V(s)  abbreviate,  respectively,  X(s,t,x,v)  and 
V(s,t,x,v).  We  immediately  have  a  uniform  bound  on  f(t,x,v), 
since 
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_d_  f(s,  X(s,t,x,v),  V(s,t,x,v))  =  0  by  Vlasov,  so 
ds 

f(t,x,v)  =  fft.Xft.t.x.v),  V(t.t.xv)) 

=  f(s,X(s,t,x,v),  V(s,t,x,v)) 

=  ffO.XfO.t.x.v),  V(0,t,x,v)) 

=  f°(X(0,t,x,v),  V(0,t,x,v)) 

=*  sup  {  f(t,x,v)  :  x  e  R,  v  e  R2  }  =  II  f°  11^ 

B.  Charge  Conservation 

Lemma  [2.1]:  JJ  f(t,x,v)dvdx  =  JJ  f°(x,v)dvdx 
Proof  :  Integrating  the  Vlasov  equation  in  v  yields 

J  3tfdv  +  J  v^xfdv  +  J  [(E+c^BMvJ'Vyf]  dv  =  0 

or  0,p  +  3xji  +  J[(E+C'1  BMv)  •  Vvf]dv  =  0 

Lemma  [2.2]:  (E+c^BMvJ'Vyf  =  Vy'^E+c^BM v)f] 

Proof:  Vv-[(E+c1BMv)f]  =  3v1[(E1+c-1Bv2)f]  +  9v2[(E2-c-1Bv1)f] 


=  (Ei+c^BvgJSv^  +  C'1fB3v1V2 

+  (E2-c-1Bvi)3v2f  -  C’1fB3v2v1 


But 


=  (E+c-1  BMv)*  Vvf  +  c-ifB(aVlv2-av2vi) 


A  A 

aVlv2  -  oy2vi  =  aVl 


( 


v2 


V 


^1+C-2(v1+V2) 


'  3v  f 


J 


Vi 


A 


V 


^1+C*2(Vi+V2) 


J 
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0  and  the  claim  is  verified. 


Using  the  compact  support  of  f  and  the  divergence  theorem  then  yields 
(2.2)  3tp  +  dxji  =  0 

Now  j!  is  compactly  supported  in  x  Vt  >  0,  as  follows  from  the  definition 
of  j  and  the  compact  support  of  f.  Integrating  (2.2)  over  x  then  gives 

3t  J  p(t,x)dx  =  0 

=>  3t  J  [  J  f(t,x,v)dv)-n(x)]dx  =  0 

=*  J  [  J  f(t,x,v)dv-n(x)]dx  =  const 

=>  J  [J  f(t,x,v)dv  -n(x)]dx  =  J  [J  f(0,x,v)dv-n(x)]dx 

=  J  [  J  f°(x,v)dv-n(x)]dx 

=>  JJ  f(t,x,v)dvdx  =  JJ  f°(x,v)dvdx  V  t  >  0,  which  proves  the 

lemma,  and  establishes  global-in-time  charge  conservation. 

C.  Field  Estimates 
(1)  Uniform  Bound  for 

We  first  establish  that  Ei(t,x)  is  uniformly  bounded  Vt  >  0  and  Vx, 
with  a  bound  that  is  independent  of  c.  Integrating  Gauss'  Law  for  Ei  with 
respect  to  x,  we  get 


and  hence 


Ei(t,x)  =  4ji  J  p(t,y)dy  +  const 


X 

Ei(0,x)  =  4n  J  p(0,y)dy 

-oo 


since  the  constant  has  been  chosen  to  be  0  already,  (the  only  choice 
resulting  in  a  finite  energy  solution).  Therefore, 

X 

E i (t,x)  =  4 n  J  p(t,y)dy 

-oo 


X 

=  4tc  J  [  Jf(t,y,v)dv-n(y)]dy 

-oo 

Using  the  non-negativity  of  f  and  lemma  [2.1],  we  have 
I  Ei(t,x)  |  <  4tcJJ  f(t,y,v)dvdy  +  4/iJ  |n(y)|  dy 

=  4kJJ  f°(y,v)dvdy  +  4tiJ  |n(y)|  dy 
=>  sup  |  Ei(t,x)  |  £  D, 

x,  t>0 

by  assumptions  on  the  data,  and  we  have  a  bound  on  E]  that  is  independent 
of  t,  x,  and  c. 

(2)  Compact  Support  of  Ei(O.x) 

Let  %  be  chosen  s.t.  f°(x,v)  =  0  and  n(x)  =  0  Vx  s.t.  |x|  >  and  consider 
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X 

Ei(0,x)=4n  J  [  Jf°(y,v)dv-n(y)]dy 

-oo 


Suppose  x  >  £.  Then 

X 

4n  J(jf°(y,v)dv-n(y))dy  =  JJ  f°(y,v)dvdy  -  J  n(y)dy  =  0 

-oo 


(2.3b) 


By  hypothesis,  the  data  terms  are  uniformly  bounded  independently  of  c, 
so  it  suffices  to  show  that 


SUP  I  J  j2('t,x±c(t-t))dx  I  is  uniformly  bounded. 

x  t  >  0  0 


.X)  =  Jv 


By  definition,  j2(t,x)  =  Jv2f(t,x,v)dv  = 


_yg_ 


v 


1  +C*2(Vi+V2) 


-f(t,x,v)dv 


j 2 ( t , x)  I  < 


jy^L 


v 


1  +C’2(v  1+V2) 


.f(t,x,v)dv 


Lemma  [2.3]:  There  exists  a  constant  D,  independent  of  c  and  depending 
only  on  the  data,  such  that 


t 


sup 

X  t  >  0 


M 


j  j 
0 


1+C’2(Vi+V2) 


f(x,x±c(t-T),v)dvdx 


<  D 


Corollary  [2.3]:  The  fields  E2  and  B  are  uniformly  bounded  :  There 
exists  a  constant  D,  independent  of  c,  such  that 

li  E2(t)  11^  +  ||  B(t)  11^  ^  D  for  all  t  >  0 
Proof  of  the  Corollary  :  By  the  integral  representations  of  E2  and  B,  we 


14 


need  to  bound 


sup 

x  t  >  o 


J  j2(x,X±C(t-x))dx  | 
0 


<  sup 

x  t  >  0 


t 

r  r 


M. 


1+C‘2(Vi+V2) 


-f(x,x±c(t-x),  v)dvdx 


<  D  by  the  lemma. 


Proof  of  the  lemma  :  We  first  derive  an  energy  identity  for  Vlasov.  The 
relativistic  energy  of  a  particle  is  c2V  1  +c_2|v|2,  so  the  energy  density  is 
given  by 


ek  -  \  c2a/1+C'2|v|2  f(t,x,v)dv 

Letting  y  =  V 1  +c2|v|2  and  differentiating  with  respect  to  t,  we  have  by 
Vlasov  and  lemma  [2.2], 

3t©k  =  J  c2  y  8tf  dv 

=  -  J  c2  y  [  v^xf  +  (E+c-1BMv)  •  Vvf  ]d  v 
=  -  9X  }  c2vifdv  - 1  c2y  [  (E+c^BmO)  *Vvf]dv 
=  -8x  J  c2v-ifdv  -  J  c2yVv-[  (E+c-iBMv)f  ]dv 
Integrating  by  parts  and  using  compact  support  of  f,  we  have 
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9tek  =  -  9X  J  c2vifdv  +  c2  J  (  [  (E+c^BMvJf  ]  •  Vvy)dv 

A 

=  -  9X  J  c2vifdv  +  c2  J  [  (E+c^BMvJf  ]  •  ^  dv 
=  -  9X  J  c2vifdv  +  J  (  [  (E+c^BMvjf  ]  •  v)dv 

Now,  since  (v  x  B)-v  =  0, 

(E+c^BMvjf  -v  =  (E+C'^v  x  B))f  *  v 
=  Ef  •  v 
=  E  •  fv 

=>  }  ( [  (E+c  ^BMvJf  ]  •  v)dv  =  J  E-fvdv 

=  E-(  J  fvdv) 

=  E-j 

so  9tek  =  -9X  J  c^fdv  +  E-j 

From  Maxwell, 

-47iji  =  9tE i  and  -4rcj2  =  9tE2  +  c9xB,  so 
j  =  -  “  (9tEi,  9tE2  +  c9xB) 
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E  *j  =  -  ^  (Ei3tEi  +  E29tE2  +  cE23xB) 

-  -£[  E'i2+ 1  E*i2)+  cE^Bl 

=  -j“|E|2  -  ^(cE23xB  +  c(3xE2)B  -  c(3xE2)B) 


(2.4) 


i  a 


i  a 


i 


=  'Sai|E|2  ‘  4^(cE2B)  +  4tcC(3xE2)B 


i 


But  again  from  Maxwell,  3XE2  =  -  ~  atB,  so  we  can  write  the  last  term  in 

c 

(2.4)  as 


^c(3xE2)B  -  -^c(13,b]b  -  -^B3,B 


■  -k- U|B|2  -  •  ^  k |B|2 


E*j  -  - 


i  k < |E|2 + |B|2)  ■  i  i  (cE2B) 


3t 


ek  = 


^fc2  |  vifdv  +  ~  E2b1  ~  k  k  <IEI2  +  IBI2) 


Regrouping  terms  gives 


(2.5) 


f-f  4nek  +  \  (|E|2  +  |B|2) 


-  S-  (4nc2  f  vifdv  +  CE2B) 


Let  e  =  47cek  +  \  (|E|2  +  |B|2)  =  4tiC2  /  V 1  +C-2|v|2  fdv  +  \  (|E|2  +  |B|2) 

and  m  =  -4jic2  J  vifdv  -  CE2B 

Then  by  (2.5), 

Be  dm 
at  ~  0x 

which  is  the  energy  identity  we  seek.  Continuing  with  the  proof  of  lemma 
[2.3],  we  integrate 


0  = 


Be 

at 


am 

ax 


over  a  backwards  characteristic  cone  with  vertex  (x,t)  and  interior  T, 
using  Green's  theorem: 


0 


JJVx.t*  (-m,e)dA 
T 


(-m,e)*nds 


x  +  ct 

=  f(-m,e)|  -(0,-1  )dy 

J  1 (y.o) 

x  -  c  t 
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t 

0 

0 

+  J  t(-m.e)|<x  x_c(t.x))-^=  (-'.OC-VTT^JjciT 

t 

which  reduces  to 

x  +  ct  t 

(2.6)  J  e(y,o)dy  =  J  [(ce-m)(x,x+c(t-x))  +  (ce+m)(x,x-c(t-T))]dt 

x  -  ct  0 

Claim  :  There  exists  a  constant  D,  depending  only  on  the  data,  such  that 
for  all  c  >1 , 

x  +  ct 

J  e(y,o)dy  <  Dc2 
x-ct 


Proof  of  the  claim  :  By  definition, 


e(y,0) 


1 

2 


<E?(y))2  +  (E°(y))2  +  (B°(y»2 


+ 


1  +C'2|v|2  f°(y,v)dv) 


^  ,  B°  ,  and  f°  are  smooth  and  compactly  supported  by  hypothesis,  while 
E,  is  compactly  supported  as  shown  earlier.  Hence  e(y  ,0)  is  compactly 
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supported  in  y  ,  and 


x+ct 


J  e(y,0)dy 

x-ct 


x+ct 

V  (E2<y»2  +  (B°(y))2]dy  + 

*/ 

x-ct 

x+ct 

4nc2  J  JV  1  +C‘2|v|2  f°(y,0))dvdy 
x-ct 


<  D  +  Dc2 

<  Dc2  for  c  >  1 


and  so  by  (2.6)  we  have 

t 

(2.7)  J  [(ce-m)(x,x+c(t-x))  +  (ce+m)(x,x-c(t-x))]dx  <  Dc2 

o 

Now  ce  ±  m  =  |  |Ei|2  (E2  +  B)2  +  47ic2J[cVl  +c*2|v|2  +  Vi  ]f(t,x,v) 


We  consider  the  integrand  in  the  last  term: 


cVl+c2|v|2  +Vi 


(c V 1  +c-2lv|2  +  Vj  )  (cV  1  +c~2|v|2  ±  vQ 

(c V 1  +C‘2|V|2  ±  Vi) 


c2  +  v2  c2  +  v| 

cVl  +C'2|v|2  ±  V!  ~  cVl  +C'2|v|2  +  | V 1 1 

9  2 

c2  +  v2 
~  2cV  1  +c-2iv|2 
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Note  that  (|v2|  -  c)2  2:  0  =>  |v2|2  -  2c|v2|  +  c2  >  0  and  so 


2  o 

v2  +  c2 
2c 


>  |v2| 


9  2 

c2  +  v2 
2cVl+C‘2|v|2 


>  tegj 

V  1+C'2|v|2 


which  implies 

(2.8)  CV1+C‘2|V|2  +  V,  >  ^=^j= 


Using  (2.6),  (2.7)  and  (2.8)  and  discarding  the  field  terms,  since  they  are 
positive,  we  have 

t  _ 

c2  J  J  [cV  1+c-2|v|2  +  Vi]  f(x,x  ±  c(t-x)),v)dvdx  <  Dc2 
oJ 


c2 


l 

JJ 


M. 


V 1  +C"2|  v]2 


f(x,x  ±  c(t-x),v)dvdx  <  Dc2 


sup 
x  t  >  o 


I 


M. 


Vi  +C’2|v]2 


-f(x,x  ±  c(t-x),v)dvdx  <  D 


for  all  x  and  for  all  t  >  0,  where  D  is  independent  of  c.  This  proves 
lemma  [2.3]. 
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Summary:  The  fields  Ei(t,x),  E2(t,x),  and  B(t,x)  are  uniformly  bounded 
for  all  x  and  for  all  t  >  0.  The  bound  is  independent  of  c  and  depends  only 
on  the  data.  These  field  bounds  enable  us  to  obtain  bounds  on  the  x- 
support  and  the  v-support  of  f. 

D.  Bounds  on  the  Support  of  f 

(1)  v  -  support:  Let  Pi  (t)  =  sup  {  |v|  :  3  x  s.t.  f(t,x,v)  *  0} 

Lemma  [2.4]:  There  exists  a  constant  D,  independent  of  c,  such  that 
PUt)  <  D(1+t). 

Proof  :  Integrating  the  characteristic  equation  for  V  gives 

o 

V(0,t,x,v)  -  v  *>  |  [E(s,X(s))  +  c-1B(s,X(s))MV(s)]ds 

t 

t 

=>  v  -  V(0,t,x,v)  +  J  [E(s,X(s))  +  C'1  B(s,X(s))MV(s)]ds 

0 

By  the  assumption  that  f°  e  cJ(R3),  there  exists  D  s.t. 

f(t,x,v)  =  f°(X(0,t,x,v),  V(0,t,x,v))  =  0  whenever  |V(0,t,x,v)|  >  D. 

A 

I VI 

Then,  since  ~~  <  1,  we  may  write 
c 

t 

(2.9)  |v|  <  D  +  f  (IIEII^  +  ||B|Ioo)ds  on  the  support  of  f(t,v). 

o 

<  D  (1+t)  by  the  uniform  field  bounds. 
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(2)  x  -  support:  Let  P2(t)  =  sup  {  |x|:  3  v  s.t.  f(t,x,v)  *  0  } 

Lemma  [2.5]:  There  exists  a  constant  D,  independent  of  c,  such  that 
P2(t)  <  D(1+t)2 


Proof:  Integrating  the  characteristic  equation  for  X  gives 


x  =  X(0,t,x,v)  + 


t 

r 


Vi(s) 


Vl+C-2|V(S)|2 


d  s 


As  in  (2.9),  using  the  compact  support  of  f°,  we  may  write 


|x|  <  D  + 


f  ^  1  +c'^ V(s) |2~ dS  0n  the  SUPPOrt  0f  f(t>''') 


s  D  +  f|V,(s)|  ds 
0 


By  lemma  [2.4],  we  have 


|x|  <  D  +  J  D ( 1  +t)ds  <  D(1+t)2  on  the  support  of  f(t,v) 
o 


E.  Bounds  on  p  and  j 

Corollary  [2.4]:  |p(t,x)|  +  |j(t,x)|  <  D(1+t)3  for  all  x  and  for  all  t  >  0, 

where  D  is  again  a  constant  independent  of  t,  x,  and  c 
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and  depending  only  on  the  data. 


Proof:  From  the  definition  of  j(t,x),  we  have  immediately: 

|j(t,x)|  <  J  | v|f(t,x,v)dv 

<  J  |v|f(t,x,v)dv 

<  J  D(1+t)dv 

v-spt  of  f 

<  D(1+t)3  for  all  x  and  for  all  t  >  0 


For  p,  we  have 

lp(t,x)|  <  J  |f(t,x,v)|dv  +  |n(x)| 

<  D(1+t)2  +  D 

<  D(1+t)2  V  x,  t  >  0, 

and  the  corollary  is  proved. 

F.  Estimates  on  Derivatives 

We  are  now  able  to  bound  the  derivatives  of  f(t),  E(t,x)  and  B(t,x), 
which  will  give  us  bounds  on  C1  norms. 

Lemma  [2.6]:  With  the  given  assumptions  on  the  Cauchy  data,  there 

exists  a  non-decreasing  function  D:  (O.oo)  (0,«>)  which  is 
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independent  of  c,  such  that 


l|f(t)llci  +  ||E(t)||ci  +  ||B(t)||c,  s  D(t) 

Proof :  We  begin  with  bounds  on  the  spatial  derivatives  of  E  and  B. 

(1)  Spatial  Derivative  of  Ej_:  Using  Maxwell,  assumptions  on  n  and  non 
negativity  of  f,  we  have 

axEi(t,x)  =  4rcp(t,x)  =  An  ( |  f(t,x,v)dv  -  n(x)) 

=>  |axEi(t,x)|  <  An  J  fdv  +  D 

<  D  -(radius  of  v-support  of  f(t,x))2  +  D 

<  D(t)  by  lemma  [2.4] 

(2)  Spatial  Derivatives  of  E?_and  B: 

To  obtain  bounds  on  3XE2  and  axB,  we  set 

K±(t,x)  =  E2(t,x)  ±  B(t,x) 

From  the  integral  representations  (2.3),  we  have 

t 

(2.10)  K^t.x)  =  K±*0(x^ct)  -  An  J j2(x,x  + c(t-x))dx 

0 

where  K±>°(x^ct)  =  E^x+ct)  ±  B°(x+ct) 
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We  estimate  3XK+  only,  since  the  bounds  for  K*  are  obtained  in  a  similar 
way.  Differentiating  (2.10)  with  respect  to  x  and  using  the  definition  of 
\2,  we  have 

t 

(2.11)  3xK+(t,x)  =  3xK+.o(x-ct)  -  4ti  J  axjv2f(T,x-c(t-t),v)dvdx 

0 

Now  |3xK+.°(x-ct)|  =  |3xE°(x-ct)  +  3xB°(x-ct)|  <  D,  since  E^,  B°  e  Cg  (R). 
As  in  [2],  we  introduce  the  following  differential  operators: 

T+:=  3t+c3x  This  is  the  derivative  along  one  of  the 
characteristics  of  the  Maxwell  system,  e.g. 

0 

T  +  f(t,X-c(t-T),V)  =  —  f(T,X-C(t-t),V) 


S  := 


3t  +  V!3x 


Note  that  by  Vlasov  and  Lemma  [2.2], 


Sf(xIx-c(t-x))v)  =  (3t  +  V13x)f(l,X-C(t-T),V) 


-((E+C‘1BMv)  •  Vvf)| 

M(x,x-c(t-t), 


V) 


=  -(VviE+c-TBMvJf) 


(T,X-C(t-x),V) 
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These  can  be  inverted  to  obtain 


cS-viT  T+-S 

3t  =  - r-  and  3X  =  — — 

C-Vi  c-v1 


Replacing  3X  in  (2.11),  we  have 


3XK+  =  3xK+.°(x-ct)  -  4n 


t 

rr  a 
Yz 


j  j 


A  (T+-S)f(x,x-c(t-x),v)dvdx 

C-Vt 


=  3xK+-°(x-ct)  -  4;t 


t 

r  r  a 

-^2— T3f 


a  (x ,x-c(t-x),v)-Sf(x,x-c(t-x),v) 
c  -  Vi  [3r 


dvdx 


=  3xK+.o(x-ct)  -  4tc 


t 

r  c  a 
V2 


3_ 

3x 


A  f(x,x-c(t-x),v)dvdx 

C-Vi 


t 


-  4ji 


r  a 

Yz 


~t-Vv((E+c-iBMv)f)| 

-V 1  '(t.x- 


dvdx 


C-V! 


C(t-T),V) 


Since  f(t,v)  is  compactly  supported,  c-vi  *  0,  and  the  v-integrals  are 
non-singular. 

r  a 
1Z. 


=>  3XK+  =  3xK+>°(x-ct)  -  4rc 


A 

C-Vi 


f(t,x,v)dv  +  4ji 


r  a 
Yz 


a  f(0,x-ct,v)dv 

C-Vi 


4  n 


t 

rr  a 
v2 


Vv-((E+c*iBMv)f)|  dvdx 

-  Vi  l(T,X-C(t-T).V) 


C-Vi 
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We  integrate  the  last  term  and  use  the  compact  v  support  of  f  to  get 


r  a 
X2 


C-Vi 


V*((E+C'1  BM  v)f)| 


(-r,x-c(t-T),v) 


dvdx  = 


t 

r  r 


J 

0 


f(E+c*1BMv)f| 

V  l(TfX-C(t-X).V)> 


(  A  \ 
YjL 


Vc-V1y 


dvdx 


So  we  have 


3xK+ 


=  3xK+>°(x-ct)  -  47i 


A 

-^-f(t,x,v)dv)  +  4k 

C-Vi 


A 

f(0,x-ct,v)dv 

C-V! 


r 

• 

(  h  \ 

+  4k 

v2 

Vv 

A 

J 

* 

lC-VlJ 

((E+C'1BMv)f)|  dvdx 


which  implies 


(2.12)  |3xK+(t,x)|  <  D  +  4k 


A 

f(t,x,v)dv  +  4k 

C-Vi 


r  a 

-^-f°(x-ct,v)dv 
C-Vi 


r 

r 

(  A  \ 

+  4k 

11x7 

v2 

||Vv 

A 

m 

lC'V1  J 

II  llE+c1  BMv| |  |f(x,x-c(t-x),v)|dvdx 


since  the  spatial  derivatives  of  E2  and  B°  are  uniformly  bounded. 
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The  v-integrals  are  taken  over  a  compact  set,  e.g.  {v:  |v|  £  D(1+t)  },  so 

A 

c-vi  *  0,  and 


(2.13)  <,  D(t)  V  t  >  0 

C-Vi 

Also,  we  have 

/  A  \  A 

1  „  A  V2  _  A 

=  - t~vw2  +  — rr^vv  i 
c-vt  (c-v^ 

and  so 


(2.14) 


/  A 

V2 


A 

C-V1 


s  -Vl|v»v2||  +  -%Tl|VvO,|| 

Ic-Vil  (c-vO 


Lemma  [2.7]:  ||  Vvvj||  s  V2  ,  i  ■  1,  2 
Proof : 


HVvOill2  . 


f  1  fv 

if  1  ^ 

Vl+c2|v|2  v( 

3J  (1+c_2|v|2)3/2 

(vjVilj 


C4(1  +C'2|v|2)' 


i  -  1  ,2;  j  -  {1  ,  2}\i 


Now 


N 


■J 


1  +  C‘2V? 


Vl+C-2|V|2  {0)  (1+c-2|v|2)3/2  (1+C'2|v|2) 


3/2 


>  0,  so 


f,  1  .  p£if 

Vl+C'2|v|2  ) 


1  Vs  <r _ 1 _ 'I2 _ 1 


(1+C2|V|2) 


3/2 


yVl  +C*2|v|2  )  1+C'2|V|2 
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so 


+ 


Then  from 


(2.15) 


llVvViH  < 


(  1 


imsli 


1+C-2|V|2  C4(1+C-2|V|2); 


<  y[2  and  the  lemma  is  proved. 


(2.13)  and  (2.14)  we  have 


I  v' 


(  A  N 
V2_ 


vc-0i. 


||  *  D(t) 


M/2 

j 


Using  (2.12),  (2.15),  and  the  bounds  on  f,  f°,  and  the  fields  then  gives 

II  dxK+(t)  <  D(t)  for  all  t  >  0,  D(t)  independent  of  x  and  c. 


A  similar  treatment  (see  remark  on  page  16  of  [2])  yields 

||3xK'(t)||oo  £  D(t),  and  the  result  is  that 
liaXK(t)||oo  <  D(t)  v  t  >  0 


Using  the  bounds  on  the  fields  themselves,  we  have 
l|E2(t)||  i  ,  1 1  B(t)  1 1  ,  <  D(t) 

C  C 

(3)  Derivatives  of  Characteristics; 

The  next  step  is  to  bound  5X,  3Vi  and  dV2,  where  d  can  be  dx,  3Vl. 
or  3v2. 

Integrating  (2.1)  gives 
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X(s,t,x,v)  =  x-J  V1(4,tfx,v)d^ 


t 

V ,  (s,t,x,v)  =  v,  -  J[E1K,X,(4))+c  'B(tX®)V2®]d^ 


I 

V2(s,t,x,v)  =  v2-  J[E2(4,X(5))-c-’B(tX(^))V,(4)]d^ 


We  first  consider  the  spatial  derivatives.  Differentiating  with  respect  to 


x,  we  find 


ax[si 

dx 


«  1  -  J 3xVi(4)d£, 


3V,(s) 

dx 


l 

=  -  J  [axE1(4,x(%)) 


c-i{B(tX(4))3x'J2ra  +  3xB(5,X(4))  ^V2(4)}]d^ 


3Vg(s) 

9x 


=  -  f  [dxE,(tX(4)) 


c-'{b (4,X(4))3xV,(y  +  :xB(|,X©)  V,($)}]dt 


Claim  :  |3 Vj(4)|  <  V2  ||8V||  where  3  can  be  3Xl  dVl,  or  3V 
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Proof : 


v 


Let  H(v)  ■  5 


Then  H  :  R2  ->  R2  and  Hj(v)  =  Vj  = 


Vi 


V 


l+C^+V*) 


i  -  1,2 


so  aVj  =  d(Hj(V))  =  vv(Hj(V))*av  =  vv(vo-av 


Ai  <  VvVill  II  av 


lavil 


(2.16) 


<  V2  1 1  aV  by  lemma  [2.7] 


|Vj| 

Taking  absolute  values,  using  (2.16),  and  the  facts  that  c>  1  and  — L  <  1, 

c 

i  =  1,2,  we  may  write 


t 

|3xX(s)|  <  1  +  J  V  2  ||3xV(5)||d$ 

s 


|3xV,(s)|  <  j{(||3xEi||c„+||dxB||J|3xX(5)|+V2||B||J|3xV(«||}d5 


t 

|3xV2(s)|  £  J {(I|3xE2|oo+  ||3xB|l_)|3xX(4)|  +V2  ||B||_  ||3xV(^)||}dE, 
s 

Using  the  bounds  on  the  fields  and  their  derivatives,  we  have 

t 

|3xX(s)|  £  1  +  J  V2  ||3xV(«||d5 

S 
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t 

|3xV,(s)|  <  J  (D(4)  |3xX($)|  +  VI  D  ||3xV«)||)d5 
s 

t 

|3xV2(s)|  S  j  (D<5)  |3«X(4)I  +  VI  D  ||3xV(4)||)d^ 

S 

We  add  these  three  equations,  and  since  se  [0,t]  and  D(4)  is  a  non  - 
decreasing  function  of  we  have  (using  the  triangle  inequality  on 

\mm ) 

|3xX(s)|  +  |3xVi(s)|  +  |3XV2(S)|  <1  + 

t 

oJ  D(t)(|dxX(4)|  +|3XV1(4)|  +|3xV2(^)|)d^ 

By  Gronwall's  lemma,  then, 

|dxX(s)|  +  |dxVi(s)|  +  iaxV2(s)|  <  etD(t)  =  D(t)  for  all  se  [0,t]. 

Now  consider  the  vi  derivatives.  Differentiating  with  respect  to  vi 
gives 

t 

3«,X(s)  -  -JdV)'/,(4)d4 

S 

t 

3v,Vt(s)  -  1  -  J[3xE,(4,X(4))3v1X(4)  +  o-'(B(4,X(4))3v1V2(4)  + 

S 

V2(4)3xB(4,X(4))3VtX(4))]d4 


t 

3v,V2(s)  =  -  J[9,E2ffi,Xffi))3v1X«)-c-'(BK,XK))3Vl01(5)  - 
s 

v1(4)axB(4,x(^)aVlx(^))]d^ 

By  the  same  arguments  used  above,  we  are  led  to 

|3v1X(s)|  +  |3v,Vi(s)|  +  |aVlV2(s)|  <  1  + 

t 

|  D(t)(|3v,X(5)|  +  |3v,V,(«|  +  |3»,V2(5)|)d5) 

0 

and  by  Gronwall, 

|3VlX(s)|  +  |3VlVi(s)|  +  |3ViV2(s)|  <  D(t),  for  all  se[0,tj 
A  similar  argument  yields  the  same  result  for  the  v2  derivatives. 

(4)  Derivatives  of  the  Vlasov  Density  f: 

The  bounds  on  the  derivatives  of  the  characteristics  enable  us  now  to 
bound  the  x  and  v  derivatives  of  f(t,v)  independently  of  c: 

3xf(t,x,v)  =  dxfoWO.t.x.v),  V(0,t,x,v)) 

=  f°(X,V)3xX(0,t,x,v)  +  VvfO(X,V)-3xV(0,t,x,v) 


34 


=»  |3xf(t,x,v)|  <  IKIL  |3xX(0,t,x,v)|  +  ||VvfO||M  ||3xV(0,t,x,v)|| 

Using  the  bounds  on  the  derivatives  of  the  characteristics  and  the 
assumptions  on  the  data  f°,  we  have 

(2.17)  ||dxf(t)IL  ^  D(t), 

Similarly, 

Vvf(t,x,v)  =  Vvfo(X(0,tlx,v),  V(0,t,x,v)) 

=>  l|Vvf(t,x,v)||  <  liaxfOlU  ||VvX(0,t,x,v)||  +  IJVvfon^  || VvV(0.t,x,v)| | 

and  so 

(2.18)  HVvfWII^  <  D(t) 

(2.17)  and  (2.18)  and  the  uniform  bound  on  f(t,x,v)  together  give 
||f(t)||c1  ^  D(t)  and  the  lemma  is  proved. 

The  t  -  derivative  of  f  is  bounded  by  using  the  Vlasov  equation,  e.g. 
atf(t,x,v)  =  -Vifx(t,x,v)  -  (E+C-1BM0)  •  Vvf(t,x,v) 

=>  |3tf(t,x,v)|  <  Ml  \U\)\L  +  (IIEWIL  +  HBWIL  ||v||)  HVvfWIL 

and  since  |vi|  <  ||v||  <  D(t)  on  the  support  of  f(t,v),  the  bounds  on  fx  and 
Vvf  yield 
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|3tf(t)|  <  D(t)  for  all  t  >  0,  with  D(t)  independent  of  c. 


Remark; 

Bounds  on  the  t  and  x  derivatives  of  j  and  p  follow  immediately  from 
the  bounds  on  the  t  and  x  derivatives  of  f  and  the  compact  v  -  support 
of  f,  using  the  definition  of  p  as  a  v  -  integral  of  f.  The  result  is 

l|3j(t)li00.  ||dp(t)||00  <  D(t)  for  all  t  >  0, 
where  3  can  be  either  3t  or  3X  and  D(t)  is  independent  of  c. 
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CHAPTER  III:  The  Modified  Problem 
A.  Modified  Field  Operators 

We  use  the  integral  representations  of  the  fields  B  and  E2: 

E2(t,x)  =  ^[E^(x-ct)  +  E^x+ct)  +  B°(x-ct)  -  B°(x+ct)]  - 

2n£®U2](t,x) 


B(t,x)  =  ^[E^(x-ct)  -  E°(x+ct)  +  B°(x-ct)  +  B°(x+ct)] 


2k e.  [j2](tiX) 


where 


£±  02](t,x)  :  =  {  [i2(x,x-c(t-x))  ±  i2(x,x+c(t-x))]dx 
0 


We  define  a  modified  field  operator,  8+  as  follows: 


(3.1) 


X  00 

e±  [g]  (t,x)  .  =  £  j g(t,y)dy  ±  ^  Jg(t,y)dy 


B.  Difference  Between  8  +  and  8+ 


By  lemma  [2.5],  there  is  a  constant  D  such  that  the  support  of  f(t,-,v)  is 
bounded  by  D(1+t)2.  Let  D0  be  any  constant  £  D. 


Lemma  [3.1]:  For  c  >  4D0  and  for  |y|  <  D0(1+t2) 
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|<e$  Dzl  -  e±  tial)|  I  s 


t  £  4  D  0/  c 

t  <  4  D  0/  c 


Proof  :  (e±  U2]  -  tC±  O2])  | 


4  1  y 

-  Jj2('t.y'C(t-x))dt  -  ~  J  j 2 ( t » rj )dr| 

0  C-oo 

t  00 

±{  Jj2(t,y+c(t-x))dT  -  -  J j2(t,T|)dr|  } 
0  y 


t  t 

=  ^(x.y-c.Ct-x^dx  ±  Jj2(x,y+c(t-x))dx 

0  0 

yoo 

i  1 

-{  ~  Jj2(t/n)dTi  ±  —  J j2(t.Tl)dTl} 

-00  y 


In  the  first  integral,  make  the  change  of  variable  ti  =  y-c(t-x)  and  in  the 
third  integral,  let  t\  =  y+c(t-x).  Then 

r  -  y 

(£±  \)z]  -  £±  [)2])|(ty)  =  J  12 ( t  -  y~~,  Tijc^dri  -  ^  Jj2(t,Ti)dTi 

y-ct  "00 


OO 

J 2  (  t  +  V~~,  Tl)(-C-1)dTl  +  ^  Jj2(t,Tl)dTl 


y+ct 
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’  y  ■  c  t 

J  [j2(l  '  y^T’  *0  ‘  i2(t,Ti)]dTi  -  ^  Jj2(t,ri)dTi 


y-ct 


y+ct 


c  f  [i2(l  +  yT“-  *0  -  j2<t,T|)]dTi  +  ^  J  j 2 ( t , rj ) dr] 

J  y  +  ct 


By  the  Mean  Value  Theorem, 


y+ct 

-  “  J 3tj2(xi,n)  ±  c  J  atj2(x2,Ti)  l^'dTi 


y-ct 


y-ct 


-  “  J  j2(t,ll)dTl  +  -  J  j2(t,Tl)dTl 

-oo  y  +  ct 


and  so  we  have 

(3.2)  |(e±[j2]  -  e±[i2] 


(t,y) 


1  y 

^  ^2  J  W2(tl,Tl)||Tl-y|dTl  + 
u  y-ct 


1  y+ct 

^2  J  |9tj2(x2,ti)|  h-y|dT]  + 


C2 


y-ct 


“  J  |j2(t,Tl)|dTi  +  ^  J  | j2 (t ,T| ) |drj 


y  +  ct 


Consider  the  first  two  integrals  on  the  RHS.  Recall,  dtj2(t,*)  is 
compactly  supported  with  support  bounded  by  a  non-decreasing  function 
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D(t).  Note  that  xi,T2  e  [0,t].  We  may  replace  the  intervals  of  integration 
with  (-D(t),  D(t))f  and  the  first  two  terms  on  the  right  in  (3.2)  are  then 


1  D(t)  1  D(t) 

£  ^  J  |3tj2(xi  ,ti)!|  Ti-y  |dTi  +  ~2  J  |3ti2(x2,ri)|n-y|dri 
°  -D(t)  °  -D(t) 


By  the  a-priori  bound  on  3tj2>  we  have 


|atj2(t,x)|  <  D(t)  where  D(t)  is  again  non-decreasing,  so  we  can 
combine  these  two  integrals  as 


1  D(t) 

<  -^  J  D(t)h-y|dTi 
c  -D(l) 


D(t) 


^  ^  J  (hi  +  lyl)dri 


-D(t) 


D(t) 

-  “c2  j  (  D(t)  +  IVi)^ 


-D(t) 


*  ^  (|y|  +  i) 


and  since  |y|  <  D(1+t2), 

,  Pill 

*  c2 

We  now  consider  the  other  2  terms  on  the  right  hand  side  of  (3.2). 
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Define  Z(t)  {x:  3  v  e  R2  s.t.  f(t,x,v)  *  0}.  Again,  by  the  a-priori 
bounds  and  the  definition  of  D0,  we  have 

(-  D0(1  +  t2)  ,  Do(1  +  t2))  3  %(l) 
where  D  and  D0  are  independent  of  t  and  c. 

Given  ye  R,  if  (y  +  cti,°°)  n  Z(  ti)  =  0  for  some  ti,  then 

(y  +  ct,  «)  n  Z  (t)  =  0  for  all  t  >  ti ,  because  the  transport  speed  of  f 

is  |vi  |  <  c.  For  the  same  reason,  if  (-oo,  y  -  ct2)  n  £(t2)  *  0  for  some  tg, 

then  (-°° ,  y  -  ct)  n  Z{\)  =  0  for  all  t  >  t2. 

4  D0 

Take  c  >  4D0  and  consider  any  (y,s)  with  — -  <  s  <  1  and  |yl  <  D0(s2  +  1). 

c 

Then  |y|  -  sc  <  |y|  -  4D0 

<  D0(s2  +  1 )  -  4 D0 

<  2D0  -  4 D0 
=  -  2D0 

<  -  D0(s2  +  1) 

So  (-oo,  |y|  -  cs)  n  Z(s)  =  0  and  because  of  transport  speed, 

(-oo,  |y|-  ct)  n  Z (t)  =  0  for  all  t  £  s.  Similarly,  (|y|  +  ct,  °°)  n  Z{ t)  =  0 

for  all  t  >  s. 

4D0 

Hence  for  any  (y,t)  with  t  >  ~~  and  |y|  <  D0(t2  +  1), 

c 

y-ct  00 

j  |j2(t,ti)|dri  »  J  |j2(t,Ti)|dn  =  0 

-oo  y  +  ct 
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and 


|(^W-e?iU)|  |  s  ^ 


4D0 

If  t  <  — -  ,  the  contribution  of  the  first  two  integrals  on  the  right  hand 
c 

D(t) 

side  of  (3.2)  is  still  bounded  by  *^r  •  There  is  now  a  contribution  from 

the  other  two  integrals,  however.  We  obtain  a  bound  on  this  by  including 
in  the  interval  of  integration  the  piece  [y-ct.y+ct]  and  use  the  bounds  on 
H2I  and  its  support  to  get 

y  - c  t  00 

~  1  Ij2(t,ri)|dri  +  “  J  |j2(t,-n)|dTi  <  |j2(t,T] ) |dr|  <  ^ 

—00  y  +  c  t 


Result: 


l<e±  Da]  -  eS  02])|(y  t)|  £ 


■m 

Uf.T*  T  "  "4V 


and  the  lemma  is  proved. 

C.  The  Modified  Problem  (RVM~) 

We  propose  a  modified  problem  in  which  the  field  operators  £+  are 
replaced  with  the  operators  £+  ,  and  in  which  the  Cauchy  data  is  the 
same: 
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(3.3) 


(RVM~) 


3tf  +  Vi3xf  +  (E+c^BMv)*  Vvf 
0  s  f«(x,v)  £  Cq(R3) 

E?(x)  &  B°(x)  e  Cq(R) 


The  data,  f°,  E«  ,  and  B°  are  considered  known.  We  also  assume  a 


smooth,  compactly  supported,  neutralizing  background  density  n(x),  as  in 
(RVM).  Charge  and  current  densities  are  given  by 


P(t.x) 


=  J  ? (t,x,v)dv  -  n(x) 


j  (t,x)  «  J  vf  (t,x,v)dv 

R2 

The  electric  and  (scalar)  magnetic  fields  are  given  by 


Ei  (t,x)  =  An 


1  p~(t, 


y)dy 


E2(t,x)  =  ^  E2°(x-ct)  +  E°(x+ct)  +  B°(x-ct)  -  B°(x+ct) 


-  2*  e+ir 2](t,x) 

B(t,x)  *  jT^(x-ct)  -  Ejtx+ct)  +  B°(x-ct)  +  B°(x+ct) 


-2k  £^r2](t,x) 
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Chapter  IV:  A-priori  Bounds  for  (RVM~) 

Assuming  existence  of  a  global,  compactly  supported  C1  solution  f ,  we 
obtain  the  following  a-priori  estimates: 

A.  Density  Estimates 


Define  X(s,t,x,v)  and  V(s,t,x,v)  by 


(4.4a) 


X(s,t,x,v)  = 


X(t,t,x,v)  =  x 


Vi(s,t,x,v) 


(4.4b) 


< 


V(s,t,x,v)  =  E(s,X(s))  +  c*1  B(s,X(s))MV1(s) 
V(t,t,x ,v)  -  v 


M  »N#  ^ 

where  X(s)  and  V(s)  abbreviate  X(s,t,x,v)  and  V(s,t,x,v).  Then 
^f(s,X(s),V(s))  =  0 


=*  f  (t.X(t).V(t))  -  Mtx.v)  »  f(0,X(0).  V(0))  =  f°(X  (0),  V(0)) 

and  since  since  f°  €  cJ(R3), 

ll'(t)IL  -  llf°IL  -  d 

(Non  negativity  of  f  (t,x,v)  also  follows  from  non  negativity  of  f°.) 
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B.  Charge  Conservation 

We  integrate  the  Vlasov  equation  in  v  to  get 


|  3tfdv  +  J  vtdxfdv  +  [(  E  +  ^p)-Vvf]d\ 

R2  R 2 


R2 


=  3tp  +  5xj  i  + 


[(ET  +  B-^)vvf]dv 


R2 


But  (E  +  c^BMv) •  Vvf  =  Vv*[F(e  +  c^BMv)],  and  since  F(t,\*) 
assumed  to  be  compactly  supported,  by  the  divergence  theorem, 


/[(  E  +  c-i  B Mv) •  V Vf]dv  =  0 

R2 


and  3tp  +  3xj  i  =  0 

We  integrate  in  x  to  get 

J  dtpdx  +  J  3xj  idx  =  J  dtpdx,  since 

J  3xj  idx  =  0  by  compact  x  support  of  j . 
=>  3t(Jpdx)  =  0  =>  Jp(t,x)dx  =  const 


J  [J  <  (t.  x,v)dv  -  n(x)]dx  =  const  =  J  [J  f  (0,x,v)dv-n(x)]dx 

=>  J  J  f  (t,x,v)dvdx  =  J  J  f°(x,v)dvdx,  which  is  global  charge 

conservation. 

C.  Field  Estimates 

1.  Uniform  Bound  on  Ei 

x  x 

Ei(t,x)  =  4 7t  J  p~(t,y)dy  =  An  J[{  f  (t,y,v)dv-n(y)]dy 

-°°  -oo 

Using  the  non-negativity  of  f  and  the  assumptions  on  the  data,  we  have 
|Ei(t,x)|  <  An  J  J  f  (t,y,v)dvdy  +  4k  J  J  |n(y)|dy 
=  An  J  J  f°(y,v)dvdy  +  4tiJ  J  |n(y)|dy 

<  D 

2.  Bounds  Qfi_E2_and_S 
First  note  that  from 

X  oo 

£±[J  231  =  “  |  j  2(t,y)dy  ±  ^  J  r2(t,y)dy 

'(t.x)  c  C  X 

lv?| 

we  get,  using  the  assumptions  on  f°  and  the  fact  that  <  1 , 

c 
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|eS  02] | (t  |  ^  l  Jir2(t,y)|dy 

^  ~  JJ  |v2|f (t,y,v)dvdy 

<  JJ  f  (t,y,v)dvdy 
=  J  J  f°(y,v)dvdy  <  D 

*N» 

It  follows  from  the  integral  representations  of  E2  and  B  and  the 
assumptions  on  the  data  functions  Eg  and  B°  (e.g.  smooth  and 

compactly  supported)  that 

|  E2(t,x)|  <  HEjIL  +  I|B°||m  +  2nD  and 

|B(t,x)|  S  ||Ej||„  +  HB°IL  +  2*D,  V  x,  t>  0 

=>  P2IL  and  ||B||  are  uniformly  bounded  in  t  and  x,  and  the 

bounds  are  independent  of  c. 

D.  Bounds  on  the  Support  of  f 

m  v  -  Support:  Let  Pi (t)  :=  sup  {|v|:  3x  s.t.  f(t,x,v)  *  0} 

Lemma  [4.1]:  There  exists  a  constant  D,  independent  of  t  and  c, 
such  that  Pi  (t)  <  D(t+1) 
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f(t,x,v)  =  ^(XfO.t.x.vJ.VfO.t.x.v))  =  0,  if  |V(0,t,x,v)|  >  D 

»v 

So  on  the  support  of  f , 

t 

|v|  <  D+  JdlEIL  +  IIBIIJds 
0 

<  D(t  +  1)  by  the  uniform  field  bounds,  and  the  lemma  is 

proved. 

(2)  x-Support:  Let  P2(t)  :  =  sup  {|x|  :  3  v  s.t.  F(t,x,v)  *  o} 

Lemma  [4.2]:  There  exists  a  constant  D,  independent  of  t,  x,  and  c,  such 
that  P2(t)  <  D(t2+1 ) 

Proof :  Integrating  the  characteristic  equation  for  X  and  using  the 

A 

fact  that  |  V  |  <  |  V  |  gives: 

t 

x  =  X(0,t,x,v)  +  J  Vi(s)ds 
o 
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t 

|X|  s  D  +  J  |v(s)|ds  on  the  support  of  f 

t 

<  D  +  J  |  V (s) |ds 


<  D  +  Dt(t+1) 


<  D(1+t2)  where  D  is  independent  of  t,  x,  and  c. 


E.  Estimates  on  Derivatives 

(1)  Spatial  Derivative  of  Ei 

X 

Ei(t,x)  =  J  p~ (t,y)dy, 

-oo 

9xEi(t,x)  =  47tp(t,x) 

»  471  (  J  f (t,x,v)dv  -  n(x)j 

|3xEi(t,x)|  <  4k  J  f  (t,x,v)dv)  +  D,  since  ne  cJ(R) 


From 
we  get 


=  4k 


J 


f  (t,x,v)dv 

|v|  <  P^t) 


+  D 


<;  4k  D*D(t+1  )2  +  D 

£  D(t+2)2,  where  D  is  independent  of  t,  x,  and  c. 
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(2)  Spatial  Derivatives  of  E?  and  B 

M  A# 

Recall  the  integral  representations  of  E2  and  B : 

E2  (t,x)  =  \  [data  terms]  -  J  j  2(t,y)dy 


1  2k 

B  (t,x)  =  2  [data  terms]  -  — 


A 

1 


-J 


j  2(t,y)dy  -  J  j  2(t,y)dy 

X 


axE2(t,x)  is  the  sum  of  x  -  derivatives  of  data  terms,  which  are 
uniformly  bounded  by  assumption,  for  example: 


|ax(BO(x-ct))[  =  iaxBO(x-ct)l  £  liaxBOH^  <  D,  since  B° «  C0(R) 
=>  l|3xE2||00  <  D,  independent  of  t,  x,  and  c 

For  axB(t,x),  in  addition  to  the  uniformly  bounded  data  terms,  we  pick  up 
an  extra  term  from 


•ax  i 


2n_ 

c 


f 


A 

1 


-J 


A 


j  2 (t , y )dy  -  J  j  2(t,y)dy 

X 


J 


4k  ~ 


j  2(t,x) 


But 


—  |j~2(t,x)|  <  71  J  |V2|f  (t,x,v)dv 


£  4k  j  f(t,x,v)dv 

v  -spt  f 
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£  4jiD(D(t+1))2  by  lemma  [4.1] 

£  D(t+1)2,  D  independent  of  t,  x,  and  c. 
Result  :  |  |3XE2|  1^  £  D  uniformly  in  t  and  x,  and 

||3xB(t)|  <  D(t+1)2  for  all  t  £  0 


(3)  Derivatives  of  characteristics 

We  obtain  bounds  on  3X,  3Vi  and  3V2  where  3  can  be  3X)  3Vl  or  3y2. 
Integrating  the  characteristic  equations  gives 

t 

X(s,t,x,v)  -  x  -  /  Vi (4,t,x,v)d4 

s 


V^s.t.x.v) 


vi 


t 


J  [E~,($,X,(5» 


A 

+  c-'B(4,X(4))V2(4)]d4 


V2(s,t,x,v) 


V2- 


J  [e'2(5.X(U)  -  c-’B(tx(«)v,(«]d4 


where  X(£)  abbreviates  X^t.x.v),  etc.  We  first  treat  the  spatial 
derivatives: 


3xX(s) 


s 


A 

dxVi<4)d$ 
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f  A 

3xv,(s)  -  -  J  {axE,fex«))axxtt))  +  c-i[B(4,x®)axv2R)  + 

S 

A 

axi(5,x®)axx(yv2(^)]}d4 

r 

axv2(s)  -  -  J{axE2(5,xK))axx(«-c-i[BU,x(y)axv1(«  + 

S 

A 

3xBft,X<y)axX(«V1K)]}dS 


A  A 


lYil  lYsl 

zs  t  rs  —  * 


Take  absolute  values  and  use  *~rL ,  l~x  £  1  (Also  note  that  s  e  [0,t]) 

c  c 


f  A 

|axx(s)|  s  i  +  J  |axv,(»|d? 

S 


|3XV,(S)|  <  J{(|axE,|  +  |3xB|)|3xX(4)|  +  c-'|B||3xV2ft)|}d4 


|3XV2(8)|  <  J{(|3xE2|  +  |3xB|)|axX(5)|  +  c-'|B||3xV1(4)|}d4 


Using  lemma  [2.7]  and  the  a-priori  bounds  on  the  fields  and  their 
derivatives,  we  may  write 


|axx(s)| 


£  1  + 


1 V  2  ||3„Vffi)||d$ 


S 
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t 

|3xV,(s)|  £  J  (D(y|dxX«)|  +  c  <DV2||SxVR)||)d5 

S 

t 

|3xV2(s)|  £  J(d(5)|3xX(5)| +  c'DV2||3,VK)||)d4 

s 

Since  we  assume  c  £  1 ,  we  may  disregard  the  c*1  term.  Also,  since  D(4) 
is  an  increasing  function  of  4.  D(4)  £  D(t)  V  ^  e  [s,t],  and  we  have 

t 

|3xX(s)|  +  |3«V,(s)|  +  |3xV2(s)|  £1  +  j  D(t)(|3xX(0l  +  l|3xV(4)||)d$ 

s 

t 

<  1  +  J  D(t)(|3xX(Ul  +  |3xVi«)|  +  |3«Vft)|)d$ 

s 

Gronwall's  lemma  then  yields 

|axX(s)|  +  |axVi(s)|  +  |axV2(s)j  <  eDWM  ^(tn-  D(t)  Vse[0,t] 

and  axX,  axVlt  and  axV2  are  all  bounded  independently  of  c. 

Now  consider  the  vi  derivatives.  Differentiating  the  characteristic 
equations  gives 

t 

3v,X(s)  =  -  J  3VlV,(4.t,x,v)d5 

s 
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and  so  by  lemma  [2.7],  we  have 


t 

(i)  |3v,X(s)|  £  J V2  ||3v,VR)||d5 

s 

t 

Then,  3v,V,(s)  -  1  -  J {3XE ,(§.XR))3»1Xtt)  +  c-'[b(5,X(4))3V)V2«) 

s 

A 

+  v2©axBfex(y)3„,X(5)]}d4 

t 

So  |3v,Vi(s)|  £  1  +  J{|3xE1R)|J3»1XR)|  +  |B(«|J3v1V2(^)|) 

s 

+  |3xB(4)L|3v,X(|)|}d^ 

|$,| 

(using  c£l  and  1  <,  1) 

c 

^  A#  0St 

Now,  using  the  a-priori  bounds  on  B,  dxB,  and  3xEi,  lemma  [2.7],  and 
the  fact  that  D(£)  is  a  non-decreasing  function  of  we  have 

t 

(ii)  |3»,V,(s)|  £  1  +  J  D(t)(|3,,X«)|  +  ||dv,V(y||)d5 

s 

Differentiating  the  equation  for  V2  gives 
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aVlv2(s) 


-  j  {d«E2ft.Xft))3v1X«)  -  c-i[B(tX(«)3„,V,(« 


V,(4)3,e(4,X(5))3v1X(5)]}d5 


and  as  before  we  may  write 

t 

(iii)  |3v,V2(s)|  S  /D(t)(|3v,X(4)|  +  ||3v,V(4)||)d4 

S 


Adding  i,  ii,  and  iii  gives 


t 

|3»,X(s)|  +  |3v,V,(s)|  +  |3v,V2(s)|  s  1  +  jD(t)(|3v,X(4)|  +  ||3v,V(4)||)d4 

s 

t 

<  1  +  J  D(t)(|3„,X(4)|  +  |3»,v,(4)l  +  |3»,V2(4)l)d4 

S 


By  Gronwall, 

|aVlx(s)|  +  jav1v1(s)|  +  |av1v2(s)|  <;  e°(t)it-si  <  eo(t)t  =  D(t)  vse  [o.tj. 
A  similar  argument  yields  the  same  result  for  the  v2  derivatives. 
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Chapter  V.  Comparison  of  Solutions 
We  have  the  original  problem,  (RVM): 

8tf  +  Vi3xf  +  (E+c^BMvJ-Vvf  =  0 
with  given  Cauchy  data 

f(0,x,v)  =  f°(x,v)f  B°(x),  and  E^x) 

all  smooth  and  of  compact  support.  The  characteristics  are  given  by 

r 

•  A 

X(s,t,x,v)  -  Vi(s,t,x,v) 

X(t,t,x,v)  =  x 

r 

V(s,tIx,v)  *  E(s,X(s))  +  c*1B(s,X(s))M  V(s) 

V(t,t,x,v)  -  v 

V 

The  modified  problem,  (RVM~)  is: 

3tf  +  v,8xf  +  (E  +  c^BMvJ-Vvf  =  0 
with  the  same  Cauchy  data 

f°(x,v),  B°(x),  and  E°(x) 

The  characteristics  for  the  modified  problem  are  given  by 
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r 


< 


X(s,t,x,v) 


Vi 


y  X(t,t,x,v)  =  x 


V(s,t,x,v)  «  E(s,X(s))  +  c*1  B(s,X(s))MV(s) 
V(t,t,x,v)  =  v 


The  solutions  of  the  two  problems  at  a  point  (t,x,v)  are 
f(tiX,v)  =  f°(X(0>t)x,v),V(0,t,x,v))  and 
f(t,x,v)  =  f°(X  (O.t.x.vJ.VCO.t.x.v)) 

Recall  that  by  the  a-priori  bounds,  there  are  c  -  independent  constants  D 
and  6  such  that  the  x-supports  of  f  and  f  are  bounded  by 
D(1  +  t2)  and  6(1  +  t2),  respectively.  Let  D0  *  max(D,  6). 

Theorem  [5.1]:  There  exists  a  non-decreasing  function 

D  :  [  0,  oo)  ->  [  0,  oo)  such  that  for  c  >  4D0  and  for  all  t  >  0, 

II  *«)  -  f  (0  IL  *  ^ 

Proof :  We  consider  the  difference  of  these  solutions  along  the 

characteristics  of  the  modified  problem.  (Note  that  if  |  X(s)  |  >  D0(1  +  s2), 
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then 


d_ 

ds 


(s.X(s).V(s)) 


(s,X(s),V(s)) 


=  0,  so  we  take  |  X(s)  |  <  D0(1  +  s2).) 
=■  ^f(s,X(s),V(s))  -  0 


=  3tf(s,X(s),V(s))  +  3xf(s,X(s),V(s))  X(s)  +  V„f(s,X(s),V(s))- V(  s) 

A 

=  3tf(s,X(s),V(s))  +  Vi(s)3xf(s,X(s),V(s))  + 


A 

[E(s,X(s»  +  cr'  B(s,X(s))MV(s)]  •  Vvf(s ,  X(s),  V(s)) 

A  A 

=■  -  {v,(s)dxf(s,X(s),V<s))  +  [e(s,X(s))  +  c  'B(s,X(s))MV(s)]  • 


A 

V„f(s,X(s),V(s))}  +  Vi(s)3*f(s,X(s),V(s))  + 


A 

[E(s,X(s»  +  c-'B(s,X(s))MV(s)]vvf(s,X(s),V(s)) 


A 


«  {[E(s,X(s))  -  E(s,X(s))]  +  c'(B(s,X(s))MV(s)  - 

A 

B(s,X(s))MV(s))}-Vvf(s,X(s,,V(s)) 

Recall  the  a-priorl  bound  I i Vvf(t,* ,*>11^  ^  D(t)  V  t£  0  where  D(t)  is 

non  decreasing.  Integrating  from  0  to  t, 
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\ 


-  <,'f)|(o.x(o,.v 


V(0)) 


J  {[E  (s,X(s))  -  E(s,X(s))  +  c-' (i(s,X(s))MV(s)  - 


B(s,X(s))M  V(s))]  •  Vvf(s ,  X(s),  V(s))}d ! 


ll(f-f)|tIL  s  Jl|Vvf(s)|| J||E(s,X(s))  -  E(s,X(s))|| 


c-’||B(s,X(s))MV(s)  -  B(s,X(s))M V(s)||}ds 


t 

S  jD(t)  {||E(s,X(s))-E(s.X(s))||  + 

0 

A  A 

c-i||B(s,X(s))MV(s)  -  B(s,X(s))MV(s)||}ds 
First  consider  ||  E(s,X(s))  -  E(s,X(s))  ||  : 

||E(s,X(s))  -  E(s,X(s))U  <  |E1(s,X(s))  -  E^s.XCs))!  + 

|E2(s,X(s))  -  E2(s,X(s))| 
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X(s) 

Now  |Ei(s,X(s))  -  E1(s,X(s))|  =  |  4n  J  (p  (s,y)  -  p(s,y))dy  | 

“OO 

X(s) 

<  4ti  J  Ip  (s,y)  -  p(s,y)|d y 

-OO 

X(s) 

<  4t:  I  J  I  f  (s.y.v)  -  f(s,y,v)|dvdy 

J  R2 

-OO 

<  An  1 1  iMs.y.v)  -  f(s,y,v)|dvdy 

<  D(t)||f(s)  -  ffs)!^  since  f(s,v)  and 
f  (s,v)  are  compactly  supported  with  support  bounded  by  a  non¬ 
decreasing  function  D(s) 

Next,  I  i2(s,X(s»  -  E2(s,X(s))  I  .  2k|  (e®[f2]-eCJ2l)|(s  -(s))| 
Triangulating,  we  write 

2re  | (e+[F 2] - e+Ua])  s  M  | ( e°[F 2] - e°D2]) |(s A\ 

*  I  (e+Lia]  -  e'fe])|(s  X)|  } 

By  lemma  [4.1],  we  have 
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For  the  other  term,  we  have 

2*  |  (  SSU  2l  -  e^Da])  |(s>S)|  *  J  I  F  2<s,y)  -  j2(s,y)  |dy 

S  ^  JJ  |v2|  |f  (s,y,v)  -  f(s,y,v)  |dvdy 

£  2k  JJ  |  f  (s,y,v)  -  f(s,y,v)  |dvdy 
<  2k  D(t)||f(s)  -  f(s)|L 

where  D(t)  is  cubic  in  t.  So  we  have  the  following  so  far: 

t 

(5.1)  ||t(t)  -  f(t)|L  s  D«>  J[D(t)||f(s)  -f(s)|L  + 

0 

A 

c-1||B(s,X(s))MV(s)  -  B(s,X(s))MV(s)| |]ds 
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and  the 


where  the  top  member  of  the  bracketed  term  applies  if  s  > 

4D0 

bottom  if  s  <  —  .  We  now  work  on  the  terms  involving  the  magnetic 
c 

fields. 


c-1| | B(s,X(s))MV(s)  -  B(s,X(s))MV(s)|| 

*  c-i||  MV(s)||  |b(s,X(s))  -  B(s,X(s))| 

|v(s)| 

Then,  since  M  is  just  a  rotation  and  <  1 ,  we  have 

c 

-  c'||S(s)|||B(s,X,<s))-B(s,X(s»| 

<  | B(s,X(s)  -  B(s,X(s))| 

Using  the  integral  representations  of  B  and  B,  since  the  data  terms 
cancel,  we  get 

|B(s,X(s))-B(s,X(s))|  =  I  (-2ne![f2]  +  2ne![j2]) |(s  x(s|)I 

-<  2"  |(ec-[F2] -£?[,,])  |(s  i(s))|  + 

2n  |(e!u2]  -  £-li2l)|(s  S(s))| 
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Again  by  lemma  [4.1]  we  have 


X  (s)  oo 

^  ~7  {  J  |  F 2(t,y)  -  j2(t,y)  I  dy  +  ||f  2 (t,y )  -  j2(t,y)  |dy  } 

X(s) 

=  “7  j  |  F 2 (t,y)  -  j2(t,y)  |d  y 
and  by  the  same  steps  as  before,  this  is 

<  2*D(t)  lifts)  -  f(s)||„ 
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A  A 

So  from  the  c1 1 1  BMV  -  BMV||  term,  we  pick  up  two  additional  terms 
under  the  integral  sign  on  the  right  hand  side  of  (5.1),  namely: 


2*  D(t)  ||f(s)  -  ^8)11..  + 


fDisl 

C2 

< 

Disl 
*  c 


Result: 


t 


** 

f 

r  d  (s)'i 
C2 

(5.2) 

l|f(t)  -  f(t)IL  S  D(t) 

J 

0 

llf(s)  -f<s)|L  +  < 

> 

D(s) 

w  c  . 

where  the  top  term  in  brackets  applies  when  s  >  4D0/c  and  the  bottom 
when  s  <  4Do/c. 


Case  1:  t<^a  .  Then  s  <  4D0/c  for  all  s  e  [o ,  t]  so  the  bottom  term  in 
c 

brackets  in  (5.2)  applies,  and 


l|f(t)  -  f(t)IL  s  D(t)  f  {  II  f(s)  -HslIL  +  ^Jds 


s  D(t)  [  J  ||f(s)  -  f(s)IL  + 


4Do/c 

r  o_m 

c 


ds  ] 
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s  £ia  + 

S  C*  + 


D(t)  J  1 1 


»(s)  -  f  (s)  ||  ds 


By  Gronwall, 


(5.3) 


1 1  f(t)  -  f(t)IL  s  ^exp(tD(t))  =  ^  for  all  t  <  ^ 


4D0  4D0 

Case  2:  t  >  — Then  from  s  =  0  to  s  =  — ,  ,  the  same  argument 
c  c 

D(t) 

applies,  and  the  error  is  bounded  by  q2  • 


For  s  e  [  4D0/c  ,  t]  ,  the  ^  error  term  aPP,ies*  and  we  have 


l|f(t)  -«(t)IL  J{  D(t)||f(s)-f(s)|L  +  ^}ds 

4Do/c 


S  ^  +  jo(t)||f(s)-f(8)ILds 


l 

-  ^  +  J  D<t)ll  f(s)-f  (S)  lloodS 


By  Gronwall,  we  again  have 


II  «(t)  -  f(t)IL  £  721 
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Result: 


for  all  t  >  0 


(5.4) 


lint)  -  t(t)iL 


m 

c2 


and  the  theorem  is  proved.  An  immediate  result  is  the  following 

Corollary  [5.1]:  Let  S{t)  =  {  x:  there  exists  v  with  either  f(t,x,v)  *  o  or 
f  (t,x,v)  *  0}.  (So  (-  D0(1  +  t2),  D0(1  +  t2))  3  S(t).)  Then 
for  all  t  >  4D0/c, 

II  E(t)  -  E(t)  ||,  +  ||  B(t)  -  B(t)||,  S  ^ 
where  ||  g(t)  ||  denotes  the  supremum  over  all  x  e  S(t)  of  g(t,x). 

The  proof  follows  immediately  from  theorem  [5.1],  the  definitions  of  E, 
B,  B,  and  the  c-independent  bounds  on  the  x  and  v  support  of  f  and  f . 
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*UD 


PART  II 


CHAPTER  VI:  Computing  Solutions 
A.  The  Program 

The  scheme  used  is  of  the  type  appearing  in  [3],  e.g.  a  particle 
method  is  used  on  the  Vlasov  equation  while  the  fields  are  advanced  by 
using  the  exact  solution  representation  applied  to  approximate  solutions. 
In  [3],  Glassey  and  Schaeffer  proved  that  such  a  scheme  for  RVM  is  first 
order  in  space  and  time.  We  take  c=1  and  begin  by  choosing  a  phase  space 
grid. 

Let  Ax,  Avi ,  AV2,  At  be  >  0. 

Define:  Ca  =  {  (x,v)  :  aiAx  ^  x  <  (ai+1)Ax,  <X2Avt  <  v-i  <  (a2+1)Av1f  CX3AV2 

V2  £  (CC3+I  )AV2} 

and  C“  =  {  (<xi+i/2)Ax,  (a2+i/2)Avi,  (a3+i/2)Av2  } 

Let  qa  =  f°(Ca)(Ax)(Avi)(Av2)  2  Jf°dvdx 

a 

C 

(  This  will  be  the  charge  of  a  particle  whose  initial  state  is  C“.  ) 

Let  Cl  =  {  a  e  Z3  :  qa  *  0  }  and  note  that  Cl  is  finite  by  hypothesis. 


Let 


e  =  Ax  and  define  S,  8E)  and  0E  by 


J1  -  |x|  if  ixl  <  1 
55  w  -  1 0  i  f  |  x  I  >  1 

5e(x)  ■  e-1S(e*1x  ) 


x 

0t(x)«  J  5e(y)dy 

-oo 


Define  a  grid  on  space-time  (  [0,®°)  x  R  )  as  follows: 

tn  =  nAt,  n  =  0, 1/2,1 ,3/2,... 
xk  =  kAx,  k  e  Z 

To  start  the  simulation,  define  approximate  quantities  (designated  for  the 
moment  by  underlines)  as  follows:  For  all  k  e  Z  and  a  e  Z3, 

E(0,xk)  =  E°(xk) 

B(0,xk)  =  B°(xk) 

X“(0)  =  (ai+1/2)Ax 
Va(t1/2  )=  V(t1/2,0,Ca) 

^  i 

To  define  the  simulation  iteratively,  assume  that  for  some  n,  E(tn,x  ), 
B(tn,xk),  X“(tn),  and  Vct(  tn+1/z)  are  known  approximations  of  E(tn,xk), 

B(tn,xk),  X(tn,0,Ca),  and  V(tn+1/2,o,Ca)  for  all  k  e  Z  and  a  e  Z3. 

At  this  point,  in  order  to  simplify  the  notation,  we  drop  the  tildas  (which 
designate  quantities  from  the  modified  problem)  and  the  underlines 
(designating  approximate  quantities).  For  the  remainder  of  this  section, 
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all  quantities  are  approximate  quantities  for  the  modified  problem. 
First,  advance  Xa  by  defining 


Xa(t)  =  Xa(tn)  +  (t-  t")V?(t"+1/2)  V  t  e  [tn,  tn+1]. 

The  next  step  is  to  compute  and  advance  the  approximate  sources  p  and  \2. 


From  f(tn+1,x,v)  =  Sqa8e(x-Xa(tn+1))5e(v-V“(tn+1/2)) 

a 

we  get  p(tn+1,x)  =  XqaSe(x-Xa(tn+1))  -  n(x) 

a 

and  j2(tn+1  ,x)  -  Eq“V“(tn+1/2)5£(x-X°(tn+1)) 

a 


Next,  the  fields  are  advanced  using  these  approximate  sources.  We  know 
p(tn+1,x)  and  j2(tn+1,x)  for  all  x.  To  advance  Ei : 


.y)dy 


-4*  J  [X  Q“8e(y-X“(tn'f'))  -  n(y)  ]dy 

.«  a 


-  4n  £qa  J  Se(y-Xa(tn+1))dy  -  4n  J  n(y)dy 

a  -  oo  -  oo 

xk 

-  4x  ^q»©e(x-X“(tn*1))  ■  4*  J  n(y)dy 

a  -oo 
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To  advance  E2  and  B  using  the  approximate  source  j2,  we  employ  the 
integral  representations  to  get 


E2(tn+\xk)  =4[E°(xk-tn+1)  +  E°(xk+tn+1)  +  B°(xk-tn+1)  -  B°(xk+tn+1)] 

-  2nj  j2<tn+1  ,y)dy 

B(tn+1  ,xk)  =^-[E2(xk-tn+1)  -  E°(xk+tn+1)  +  B°(xk-tn+1)  +  B°(xk+tn+1)] 

xk  00 

-2n[  /j2(tn*1,y)dy  -  J  j2 (tn+ 1  ,y)dy  ] 

-°°  xk 

The  data  terms  are  known.  To  treat  the  integrals,  we  proceed  as  follows: 
2rcJ  j2(tn+1  ,y)dy  =  2it  J[£q«V“(t"*,'2)8e(y-X“(tn*,))]dy 

a 

-  2k][>J[  V“(tn+,/2)8t(y-X“(tn*'))]dy 

a 

■=  2Ii5>V“(t"t,,2)J’  8t(y-X“(tl’*’))dy 

a 

=  2n£q“V“(tn+1/2) 
a 

and  for  the  integrals  in  the  expression  for  B,  we  have 
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and 


x 


k 


Jj2(tn+1,y)dy 

-oo 


oo 

Jj2(tn+\y)dy 

k 


-  J[S<l“V“(t,'+1/2)8e(y-X“(tn*1 ))] d  y 

-  oo  CX 


*  [Xq“V2(tn*1,2)]©s(x-xa(t"*1)) 

a 


oo 

-  /[2 q“V“(tn+1/2)Se{y-X“(tn*1  ))]d  y 
k  a 


oo 

-  £q“V“(t"*,/2)  J 8t(y-X“(tn*'))dy 

“  xk 

-  5>v“(t"*,,2)[i  -  0£(x-x“(tnt1))] 

a 


The  result  is  that  the  fields  are  advanced  using 


Ei(t"*\xk)  =4n£q«0e(x-X“(r’))  -4k  J  n(y)dy 

a  -oo 


E2(tn+\xk)  =HE2(xk-tn+1)  +  E°(xk+t"+1)  +  B°(xk-tn+1)  -  B°(xk+tn+1)] 

-  2K£q«v“(tn+1/2) 

a 

B(tn+1,xk)  =^-[E2(xk-tn+1)  -  E°(xk+tn+1)  +  B°(xk-tn+1)  +  B°(xk+tn+1)] 

+  27rXq«V“(tn+1/2)[l  -  20£(x-Xa(tn+1 ))] 
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We  now  define  the  fields  E(tn+\x)  and  B(tn+1,x)  V  x  e  R  by  linear 
interpolation: 

B(tn+1,x)  =  X  B(tn+1  ,xk)eS  (x-xk) 
k  e 

with  similar  expressions  for  Ei  and  E2.  The  final  step  is  to  advance  the 
momenta.  We  take 

V“(tn+3/2)  =  V“(tn+1/2)  +  AtE1(tn+1,Xa(tn+1))  + 

yB(tn*\x0(tn*,))[v£(tn+’/2)  +  V“(tn*3'2)] 

Similarly, 

V“(tn+3'2)  -  v“(tn*''2)  +  AtE2(tn*1,X°(tn*1))  - 

yB(t"+11X°(t"*,))[v?(t 

Or,  as  a  system, 

Va(tn+3/2)  =  Va(tn+1/2)  +  AtE(tn+1,Xa(tn+1))  + 

yB(tn*^Xa(tn*,))M[v“(ln<,'2)  +  V“(tn*3'2)] 

This  nonlinear  system  is  always  uniquely  solvable  for  Va(tn+3/2)  and  this 
completes  one  step  of  the  scheme. 

The  program  was  tested  extensively,  first  by  comparing  against  single 
particle  problems  for  which  exact  solutions  could  be  computed,  and  then 


n+1/2)  +  V?(tn+3/2)] 
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by  generating  steady-state  solutions  (  see  appendix  A  )  which  also 
provided  "exact"  solutions  against  which  the  program's  output  could  be 
directly  compared.  We  found  the  scheme  to  be  better  than  first  order 
accurate  in  x,v,  and  t,  but  not  quite  second  order. 
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CHAPTER  VII:  Numerical  Experiments 

A.  Scaling 

The  analysis  of  the  convergence  of  solutions  of  (RVM)  and  (RVM~)  as 
a  function  of  the  (increasing)  parameter  c  was  necessarily  done  with  c 
appearing  explicitly.  We  hoped  to  demonstrate  this  convergence  by  direct 
comparison  of  computed  solutions  of  the  particle  scheme  coded  for  the 
modified  problem  with  those  produced  by  the  Glassey/Schaeffer  scheme 
from  [3],  given  the  same  Cauchy  data.  A  complication  is  that  in  [3],  the 
speed  of  light  is  taken  to  be  1  (as  are  the  charge  q  and  the  electron  rest 
mass  m)  and  in  order  to  compare  the  output  of  the  two  programs  when 
provided  with  the  same  data,  we  re-scale  the  modified  problem  to  achieve 
c  =  1 .  (See  appendix  B)  The  analog  of  the  parameter  c  becoming  large  in 
the  unsealed  problem  then  is  the  v-support  of  the  data  (and  solutions) 
becoming  small  in  the  scaled  problems. 

B.  Data 

After  preliminary  tests  on  simplified  problems  (for  instance  problems 
with  E2  =  B°  =  0),  we  analyzed  results  with  2  basic  sets  of  Cauchy  data, 

providing  a  "symmetric"  problem  and  an  "asymmetric"  one. 

1.  Symmetric  Problem 

For  the  first  problem,  the  data  is  as  follows:  Define 

(7.1)  f°(x,v1,v2)  =  [(1-x2)(1-v12)(1-v22)]3  for  |x|  <  1,  |vt |  <  1 ,  |v2|  <  1 

and  f°  =  0  otherwise. 
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% 


Then,  since 


1 

J(1-S2)=>dS»f§. 

-1 


we  see  that 


J  f°dv  =  (§§J(1-x2)3  for  jx|  <  1  and  0  for  |x|  £  1 
Define  n(x)  -  ^|J(1-x2)3  +  2x(1-x2)5 


Then  the  condition  of  global  neutrality  is  met,  since  the  first  term 
in  n  cancels  with  J  ,f0dv  and  the  second  term  is  odd,  hence  has 

integral  zero.  Note  also  that  the  smoothness  requirements  on  the 
data  are  met. 

As  data  for  E2  and  B,  we  choose 

4  -  a,(1-x2)3 
B°  .  a2<1-x2)3 

where  the  parameters  ai  and  a2  allow  us  to  adjust  the  size  of 
the  data. 

The  above  is  for  the  unsealed  problem,  i.e.,  the  case  in  which 
the  data  are  fixed  and  c  is  allowed  to  increase.  For  purposes  of 
computation  and  comparison,  we  scale  to  c=1  and  allow  the 
support  of  the  data  to  decrease.  As  discussed  in  appendix  B,  we 
take,  using  bars  to  designate  quantities  in  the  c=1  problem, 
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fo  (x,v)  =  c2f°(cx,cv) 


-  c2[(1-c2x2)(1-e2v12)(i-c2v22)]3  for  JxJ,  |v-i|,  |v2|  <  ^ 
and  f°  =  0  otherwise 

Eo2  (x)  -  f  (1-C2X2)3 


B°(x)-f  (1-c2X2)3 

Also, 

n(x)  =  n(cx) 

=  (§§J  (1-c2x2)3  +  2cx(1-c2x2)5  if  |x|  <  ~  and 
n(x)  =  0  otherwise 

It  is  clear  that  allowing  c  to  increase  in  the  unsealed  problem 

corresponds  to  the  support  of  the  data  in  the  scaled  problem 

decreasing  like  ~ . 

c 

2.  Asymmetric  Problem 

In  the  second  problem,  the  data  functions  above  are  shifted  as 
follows: 

We  translate  f°  as  given  in  (7.1)  by  1  in  all  three  coordinate 
directions  in  phase  space. 
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f°(x,v)  -  [(1-(x-1)2)(1-(v,-1)2)(1-(v2-1)2)]3 
if  0  <  x,  Vi,  v2  ^  2  and 
f°  =  0  otherwise 

2 

Then  J  [l -(s-1  )2]3ds  =  ||  and  we  define 
o 

n(x)  -  r||f(1-(x-1)2)3  +  2(x-1)(1-(x-1)2)5 

and,  as  in  the  symmetric  problem,  we  have  global  neutrality. 

In  a  similar  way,  we  set 

E*  -  a,(1-(x-1)2)3 
B°  =  a2(1-(x-1)2)3 

In  this  case,  the  initial  conditions  provided  for  the  scaled  problems 
are 

r°(X.V)  =  C2[(1-(CX-1)2)(1-(CV,-1)2(1-(CV2-1)2)]3 

2 

if  0  <  x,  Vi,  v2  <  ~  and  0  otherwise. 

o 

Eo  (x)  -  |1[1-(cx-1)2]3 

B°  =  ^[l-(cx-l)2]3 
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Also, 


n  (x)  =  pj  [l-(cx-1)2]3  +  2(cx-1)[l-(cx-1)2]5 


if  0  <  x  < 


2 

c 


and  n  (x)  =  0  otherwise. 
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Chapter  VIII:  Summary  of  Results 
A.  Expected  convergence  rates 

In  attempting  to  observe  the  convergence  rate  of  the  solutions  as  c 
increases,  we  compared  the  values  of  the  fields  E2  and  E2  or.  the  union 

»N# 

of  the  support  of  f  and  the  support  of  f ,  i.e.  where  there  is  charge. 
(Recall  that  E2(t,*)  is  not  compactly  supported.) 

We  chose  values  of  c  of  25,  50,  100,  and  200,  and  compared  the 
outputs  of  the  two  programs  at  t  =  0.08.  Since  the  scheme  from  [3] 
operates  "at  the  CFL  boundary",  e.g.  Ax  =  At,  a  balance  had  to  be  struck 
between  a  small  enough  Ax  (which  determines  the  number  of  particles) 
for  good  resolution,  and  a  reasonable  number  of  timesteps  to  avoid 
excessive  run  times.  We  settled  on  total  particle  numbers  of  64,000  for 
the  modified  program  and  59,319  for  the  unmodified  one,  (The 
discrepancy  is  due  to  a  slight  difference  in  the  way  the  programs 
initialize  the  particles)  with  corresponding  Ax's  ranging  from  0.002  for 
the  c=25  runs  to  0.00025  for  c=200.  With  At  =  Ax,  this  resulted  in  10 
timesteps  for  c=25,  increasing  to  320  timesteps  for  c=200. 

We  were  interested  in  documenting  the  following  3  main  results: 

1 .  Convergence  of  solutions  of  the  two  problems  as  c  grows. 
According  to  corollary  [5.1],  we  should  observe 
||  E2(0.08)  -  E2(0.08)  ||s  <  ^ 

This  result  applies  to  the  unsealed  problem.  For  the  scaled 
problems,  we  have 

E2(t,x)  =  c-1E2(t,cx) 

E2(t,x)  =  c*1  E  2(t,cx) 
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so  II  E2(t)  -  E2(t)  Hs  =  ||  cE  2<t)  -  cE2(t)||s 
II  E2(t)  -  E2(t)||,  s  ^ 

2.  Avoidance  of  the  CFL  restriction  in  the  modified  problem, 

e.g.  attainment  of  comparable  agreement  with  the 
program  for  the  unmodified  problem  when  the  modified 
program  uses  much  larger  timesteps. 

3.  Improvements  in  accuracy  in  solutions  of  the  modified 

problem  over  those  computed  without  the  data  terms 
present  in  the  solution  representation.  A  significant 
difference  here  would  demonstrate  the  value  of  the 
appearance  of  the  data  terms  involving  the  original  Cauchy 
data,  which  again  would  not  be  present  were  the  Darwin 
and  quasi-electrostatic  modifications  of  Maxwell's 
equations  made  first. 
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B.  Results 


We  first  give  results  for  the  symmetric  problem. 


c 

At 

max  IE2-E2I 

*«<• 

with  At  =At 

max  IE2-E2I 

with  large  At 

max  IE2-E2I 

with  no 

data  terms 

25 

0.002 

3.514E-03 

3.532E-03 

1 .347E-02 

3.8 

50 

0.001 

7.120E-05 

49.3 

7.295E-05 

1.680E-03 

8.0 

23.6 

1  00 

0.0005 

9.057E-06 

9.471  E-06 

2.102E-04 

8.0 

23.2 

200 

0.00025 

1.135E-06 

8.0 

At  =  0.004: 
1.180E-06 

At  =  0.008: 
9.546E-06 

2.629E-05 

8.0 

23.2 

Remarks  :  1 .  Convergence  rate  -  Column  4  shows  that  as  c  is  doubled,  the 


solutions  converge  at  the  rate  of  1/c2,  as  predicted  in 
Chapter  V. 

2.  Avoidance  of  CFL  limitation  -  The  value  of  At  used  in  the 

c=25,50,100  runs  was  0.008,  which  in  the  modified  scheme 
yielded  the  same  accuracy  as  the  CFL-limited  unmodified 
scheme.  For  the  c=100  runs,  this  gives  a  factor  of  16  in  the 
size  of  the  timestep.  In  the  c=200  run,  we  achieved  the 
same  accuracy  with  At  =0.004,  which  is  again  a  factor  of 
16  better  than  the  unmodified  scheme. 

3.  Inclusion  of  data  terms  -  Column  8  shows  a  factor  of 

approximately  23  in  the  accuracy  of  the  solutions  when  the 
data  terms  are  present  in  the  solution  representation  vs 
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when  they  are  omitted. 

The  analagous  results  for  the  asymmetric  problem  are 


c 

At 

max  IE2-E2I 

with  At=At 

max  IE2-E2I 

with  large  At 

max  IE2-E2I 

with  no 

data  terms 

col6 
co  13 

25 

0.002 

4.603E-03 

4.630E-03 

1.443E-02 

3.1 

50 

0.001 

1.337E-04 

33.4 

BBSl 

1.762E-03 

8.2 

12.8 

100 

0.0005 

1 .588E-05 

8.7 

1.650E-05 

2.177E-04 

8.1 

13.7 

200 

0.00025 

1.916E-06 

8.3 

At  =  0.004: 
1.972E-06 

At  =  0.008: 
1.016E-05 

2.706E-05 

8.0 

14.1 

Remarks  :  1.  The  1/c2  convergence  rate  is  again  shown  by  the  entries  in 


column  4. 

2.  Allowable  At  is  again  16  times  as  large  as  in  the 

unmodified  problem. 

3.  Although  not  as  great  as  for  the  symmetric  problem,  the 

increase  in  accuracy  resulting  from  the  inclusion  of  the 
data  terms  is  still  better  than  an  order  of  magnitude. 
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Appendix  A:  Steady  State  Solutions  for  (RVM) 

We  seek  a  solution  of  the  time  -  independent  Vlasov  equation  (with  c  =  1) 

(A.l )  vi3xf(x,v)  +  (E(x)  +  c*1B(x)Mv)-Vvf(x,v)  =  0 

We  impose  the  condition  E2  =  0.  Since  in  the  time  -  independent  (and  1 
space,  2  momenta)  case,  Faraday's  Law  of  Induction  becomes 

3xE2  =  0 

=»  E2  =  const,  and  E2  ■  0  is  the  only  finite  energy  solution.  The 

equations  of  the  characteristics  are  then 

•  A 

X  -  V, 

V|  =  Ei(X)  +  V2B(X) 

V2  =  -VtB(X) 

Since  Et  is  an  electrostatic  field,  there  is  a  potential,  U,  such  that 
E  i  (x)  =  -U(x) 


We  introduce  (B  ,  an  anti  -  derivative  of  B,  e.g. 


Consider  the  quantities 


£  :=  Vl+|v|2  +  U 
t  :=  V2  +  (B 


Along  characteristics  of  the  time  -  independent  Vlasov  equation,  these 
quantities  are  conserved: 


^S(X(t),V(t))  =  dx£x+vv£-v 


=  u'(X)X  + 


v-v 

Vl+|V|2 


=  Vi-Vt  +  V2*V2  +  U'(X)X 


=  Vi[  -u'(X)  +  v20B  (X)]  +  v2(-v1(B  (X)) +  u,(X)v1 
=  0 

^(v2(t)  +  0(X(t))  =  v2  +  GB'(X)X 

»  -v:( B(X)  +  (B  (xjv, 

=  o 

It  follows  that  f(x,v)  :■  g(S,£)  is  a  solution  of  (1)  for  any  (sufficiently 
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smooth)  function  g. 


The  potential  U  satisfies  Poisson's  equation: 

-  U  =  4np  =  4ti(  Jfdv  -  n)  or 

(A. 2)  -  U  (x)  -  An  [  J  g (V 1  +|v|2  +  U(x),  V2  +  GB  (x))d v  -  n  ( x ) J 

In  the  time  -  independent  case,  Ampere's  Law  becomes 
-3XB  =  4nj2,  so  GB  satisfies 

-  GB  =  4rcj2  ■  4k  J  V2fdv,  or 

(A. 3)  -  GB  (x)  =  4tc  J  V2g(Vl+|v|2  +  U(x),  v2  +  GB  (x))dv 

We  have  then  equations  satisfied  by  U  and  GB  and  the  task  is  to  find 
U,  GB  ,  g,  and  n  so  that  (A.2)  and  (A. 3)  are  satisfied,  neutrality  holds, 

and  the  appropriate  requirements  for  compact  support  and  (to  the  extent 
possible)  smoothness  are  met. 

For  simplification,  we  will  require  (A.2)  and  (A.3)  to  hold  on  (0,1).  The 
various  functions  appearing  will  be  defined  on  (-1,0)  as  even  or  odd 
extensions,  and  outside  (-1,1)  as  constants.  With  these  requirements  in 
mind,  we  impose 

(BC)  U'(0)  =  U'(1)  =  GB  '(0)  -  GB  '(1)  =  0 
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U  and  (B  are  extended  evenly  on  (-1,0)  and  as  constants  on  R\(-1,1) 


as  follows: 


U(x)  =  U(1)  =  0,  |x|*  1 

U(-x)  =  U(x) 

(B(x)  =  0(1),  |x|  a  1 
(B  (-x)  -  (B  (x) 

For  n,  we  require 

n(x)  =  0,  |x|  £  1 

n(-x)  =  n(x) 

On  g,  we  impose  the  condition 
g(S,i)  =  0  if  S  >1 

Note  that  this  ensures  the  compact  x-support  of 
|x|  >  1  =>  U(x)  =  0 

=>  £  =  Vl+|v|2  +  U  £  1 

=>  g(£,£)  f(x,v)  *=  0 


f,  since 
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With  these  definitions,  (A.2)  and  (A. 3)  are  satisfied  V  x.  Note  also  that 
the  neutrality  condition  is  built  into  the  boundary  conditions,  since 


J  pdx  -  1/4*  f  E'dx  =  ^  (E,(1)  -  E,(-1» 

-  ^(u'(-i)-u'd))  =0 


We  now  proceed  to  choose  g.  Let 

s(S.fc)  -«-2I(_.i)<£) 

=  (v*  +  2v2(B  +  IB 2)  i(._  ,,(£) 

With  this  choice,  we  find 

J  V2fdv  =  |  v2(v2  +  2v2(B  +  (B  2)l( _  1)(S)dv 

Now  £>  =  Vl+|v|2  +  U  <  1  <=>  Vl+|v|2  <  1  -  U 

<*  1+|v|2  <  (1-U)2 
|v|2  <  (1-U)2  -  1 
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<=>  |v|  <  V(1-U)2  -  1 


Note  that  if  U  ^  0,  this  is  always  defined.  The  integral  becomes 


V2(V2  +  2V2(B  +  (B  )dv 

jv |  <  *J(1-U)2-i 


J  Vi+M2 

|v|<^(i-U)2-i 


=  2tc(B  |  (1-U)3  -  (1-U) 


and  the  equation  for  (B  becomes 


(A-4)  -  ~  ©  (X)=^[((1-U(X))3  .  3(1  -U (x))  +  2]©  (X), 


when  U(x)  <  0. 


We  will  also  need  to  evaluate  J  fdv: 


J  fdv  =  J  (V*  +  2v2(B  +  (B  2)I(_  (£)dv 


J  (v2+(B2)dv 

|v|  <  ^  (i-U)2- 1 
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which  gives 


(A. 5) 


Jfdv  =  |  ((1-U)2-l)2  +  2((1-U)2-i) 


Next,  we  choose 


U(x) 


-  A  (  1  -  X  2  )  3  if  |  X  |  <  1 
.  0  if  j  x  |  >  1 


where  A  is  a  parameter  whose  purpose  will  be  explained  shortly. 
Inserting  this  into  (A.4),  we  have 
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(A. 6)  Go  (x)  =  -  ~Y~  g(x)Go  (x)  where 


a(x)  =  1 


(l+A(1-X2)3)3 
0  if  |  x  |  > 


3(l  +A(1  -x2)3)  +  2  if  |x|  <  1 


Comment:  Different  choices  of  g  can  simplify  the  equation  for  .  For 
example,  if  g  =  (C+tS»)I(  _  (£)  where  C  is  a  constant  chosen  to  force 

-  99 

non-negativity  of  f,  we  obtain  (using  the  same  U),  Go  (x)  =  I(x)  where 

Z  is  a  16th  degree  polynomial  in  x.  This  can  be  integrated  twice  to  find  a 
closed  form  solution  £B  (x).  It  turned  out,  however,  that  this  solution 

resulted  in  excessive  amounts  of  charge  present  to  maintain  the  steady 
state.  Computation  of  the  fields  involves  subtracting  an  integral  of  the 
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background  ion  density  from  an  integral  of  the  electron  density.  Since 
these  turned  out  to  be  large  numbers,  an  extremely  fine  refinement  of  the 
mesh  (resulting  in  a  prohibitively  large  number  of  particles)  was  required 
to  resolve  the  fields.  The  above  choice  of  g  was  deemed  the  best  after 
much  experimentation.  To  compute  a  solution  of  (A.6),  we  use  a  modified 
shooting  method.  We  choose  (arbitrarily)  OB  (1)  =1  and  require  GBO)  =  0. 

We  solve  an  initial  value  problem  with  this  data  specified,  and  vary  the 
parameter  A  until  we  achieve  a  zero  slope  at  x  =  0.  GB  is  then  extended 

as  described  at  the  beginning  of  this  section. 


The  difference  scheme  used  is 


(B  kk1  -  2©  k  +  (3 


k-1 


87l2 


(Ax)2 


ok(x)GB 


or  OB  k'1  =  2GB  k  -  OB  k+1  -  (Ax)2  o(xk)  OB 


where  GB  n+1  =  OB  n  =  1,  Ax=^,  and  we  try  to  achieve  OB  _1  =  OB1 

With  n  =  500,  the  value  of  A  required  to  achieve  the  boundary  conditions 
is  found  to  be  A  =  0.808195.  At  this  point,  we  consider  OB  a  known 

function.  In  practice,  we  have  OB  at  500  mesh  points  in  the  interval 
[0,1].  In  the  program,  where  required,  values  of  OB  at  intermediate 
points  are  linearly  interpolated,  while  values  of  B  =  GB  '  are  obtained  using 
the  Mean  Value  Theorem. 

Having  found  OB  ,  we  are  able  to  determine  n.  For  |x(  <  1,  we  have  from 
(A. 2): 
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1 


n(x)  =  J  g(Vl+|v|2  +  U,  V2  +  OB  )dv  +  —  U  (x) 

ff 

Using  (A.5)  and  substituting  in  for  U  and  U  ,  we  have 

n(x)  =  |  [(l+A(1-x2)3)2  -i]2  +  x(B2(x)[(i+A(1-x2)3)2-i] 

-  ^  (1-x2)(Sx2-1)  for  |x|  <  1 
and  n(x)  =  0  for  |x|  >  1 

Note  also  that  the  finite  energy  solution  of  the  steady  state  RVM  problem 
also  solves  the  steady  state  version  of  the  modified  problem,  (RVM~), 
with  the  condition  E2  =  0.  To  show  this,  we  write  steady  state  (RVM~) 
in  integral  form: 

W  +  (E  +  BMv)*  V  vf  =  0 
x 

Ei(x)  =  4 n  J  p~ (y)dy 

-00 

E2(x)  =  -2nj  r2(y)dy 

X  00 

B(x)  =  -2k [  J  f2(y)dy  -  J  J~2(y)dy] 
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The  requirement  E2  ■  0  means 


Appendix  B:  Rescaling  to  Achieve  c=1 

Suppose  (f,E,B)  is  a  solution  of  (RVM)  for  some  value  of  c,  i.e. 
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dtf(t,x,v)  +  Vi  3xf(t,x,v)  +  (Ei(t,x)  +  c*1  v2B(t,x))aVlf(t,x,v) 

+  (E2(t,x)  -  c-1v1B(t,x))3v2f(t,x,v)  =  0 
where  p(t,x)  =  J  f(t,x,v)dv  -  n(x) 

j(t,x)  =  J  vf(t,x,v)dv 

A  V 

V  Vl+C-2|V|2 
x 

E,(t,x)  »  4k  j  p(t,y)dy 

-  oo 

atE2(t,x)  =  -caxB(t,x)  -  47ij2(t,x) 
atB  =  -caxE2(t,x) 

We  re-scale  as  follows: 

Let  t  =  t,  x  =  cx,  v  =  cv  and  define: 

f(t,x,v)  =  c2f(t,x,v)  =  c2f(t,cx,cv) 
n(x)  =  n(x)  =  n(cx) 
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p(t,x)  *  J  f  (t,x,v)dv  -  n(x) 

=  J  c2f(t,cx,cv)dv  -  n(cx) 

=  J  c2f(t,cx,w)c-2dw  -  n(cx) 

-  J  f(t,cx,w)dw  -  n(cx)  =  p(t,cx) 

A  A  V 

j  (t,x)  =  J  V  f  (t,x,v)dv  where  0  :  =  ^== 

A 

=  J  v  c2f(t,cx,cv)dv 


Note  that 


so 


=  |  vc2f(t,cx,w)c-2d  w 

.  .A  CV  A 

<cv)  °v' 

s 

v  =  c'Mcv)  ,  and  since  w  =  cv, 

A 

v  =  C'1  w  and  we  have 
j(t,x)  =  J  c-1wf(t,cx,w)dw 
=  c1  j  (t,cx) 
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Now  define 


x 

Ei(t,x)  =  4k  J  p(t,y)dy 

-oo 

X 

-  4k  J  p(t,cy)dy 

-oo 

CX 

=  4n  J  p(t,z)c-1d  z 

-oo 

CX 

=  c -Hk  J  p(t,z)dz 

-oo 

=  c*1E1(t,cx) 

Similarly,  define 

E2(t,x)  =  C‘1E2(t,cx) 

B(t,x)  =  c1B(t,cx) 

Then  (  f  ,E,  B)  solves  (RVM)  with  c  =  1.  To  show  this,  we  first 

demonstrate  that  E2  and  B  are  solutions  of  the  one  dimensional 
Maxwell  System: 
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atE2(t,x)  +  axB(i,x)  =  c-1[atE2(tfcx)  +  caxB(t,cx)] 

=  c1  [-4nj2(t,cx)] 

«  -4jic-1[cj2(t,x)] 

-  -4tc  j  2(t,x) 

atB(t,x)  +  axE2(t,x)  =  c*1  [atB(t,cx)  +  caxE2(t,cx)] 

=  0 


We  also  have  axEi(t,x)  =  4rcp(t,x)  by  definition,  and  it  remains  to  show 
that  the  Vlasov  equation  with  c  =  1  is  satisfied.  Toward  this  end,  we 
first  compute 

at  f(t,x,v)  =  3,[  c2f(t,cx,cv)]  =  c2atf(t,cx,cv) 
axf(t,X,v)  =  ax[c2f(t,CX,CV)3  =  c3axf(t,CX,CV) 
av.f(t,x,v)  =  aVi[c2f(t,cx,cv)]  =  c3aVjf(t,cx,cv)I  i  =  1,  2 
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=>  3tf(t,x,v)  +  Vi3x  f  (t,x,v)  +  (E1(t,x)+v2B(t,x))8v1  f  (t.x.v) 

A 

+  (E2(t,x)-viB(tIx))9v2  f  (t.x.v) 

=  c23tf(t,cx,cv)  +  c-^cvi)  c33xf(t,cx,cv) 

+  (c*1  Ei  (t,cx)  +  c*1(cv2)A  c_1  B(t,cx))c33v1f(tIcx>cv) 
+  (c1  E2(t,cx)  -  c-^cvi)*  c-1B(t,cx))c3aV2f(t>cx,cv) 

=  c2atf(t,CX,W)  +  C2w3xf(t,CX,W) 

+  (c^E^t.cx)  +  c*2w2B(t,cx))c3aWlf(t,cx,w) 

+  (c1  E2(t,cx)  -  c-2wiB(t,cx))c3aW2f(t,cx,w) 

=  c2{(3tf  +  w3xf  +  l(E  +  c-iBMw)*Vwf)]}  |(t  cx  w) 

=  0 

Finally,  we  compute 

3tEiM  -  3,[c  'E,(t,cx)] 

=  c-^tE^t.cx) 
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=  c-1[-4nji(t,cx)] 


=  ~4tc[c*^  ji  (t,cx)J 


=  -47t’ji(t,x),  and  as  claimed, 

=  (  f ,  E,  B)  solves  (RVM)  with  c  =  1 . 

We  note  here  that  the  modified  problem  admits  the  same  scaling.  To  see 
this,  suppose  (F,  E,  B  )  solves  (RVM'-)  for  some  value  of  c  >  1.  Then  the 
Vlasov  equation  and  the  expressions  for  p  and  j  are  the  same  as  those 
in  (RVM)  with  f,  E  and  B  replaced  by  F,  E  and  B.  In  (RVM~),  the 
fields  are  given  by 


E2(t,x) 


B(t,x) 


^[E°(x-ct)  +  E^x+ct)  +  B°(x-ct)  -  B°(x+ct)] 

2  n  f  ~ 

VJ  i  z(t,y)dy 


^-[E^x-ct)  -  Eg(x+ct)  +  B°(x-ct)  +  B°(x+ct)] 


a 

2k  f  2k  f  ~ 

~  J  j  2(t,y)dy  +  ~  J  j  2(t,y)dy 


We  define  the  following: 
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f(t,x,v)  :=  C2f(t,CX,CV) 


n(x)  :  =  n(cx) 

P(t,x)  :=  p  (t,CX) 

r  (t|X)  :=  C-l[(t,CX) 

E(t,x)  :=  C'^  E  (t,cx) 

B(t,x)  :=  C'1  B(t,CX) 

With  these,  we  find 

fe(t.x)  =  C'1  E  2(t,cx) 

=  [EgCcx-ct)  +  E^cx+ct)  +  B°(cx-ct)  -  B°(cx+ct)]  - 

2ti 

~2  J  j  2(t,y)dy 

Since  Eg(x)  =  ^(O.x)  =  c-’E^O.cx)  =  c_1eJ(cx),  and  a  similar  result 
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holds  for  B(x),  we  have 


E2(t,x)  -  |[c-'E°(c(x-t))  +  c  'E°(c(x+t))  +  c-'B°(o(x-t))  -  o  ' B°(c(x+t))] 

2it  f  - 

-  —  J  c-'j2(t,y)dy 

=  |[Ej(x-t)  +  l°2(x+t)  +  BP(x-t)  -  BP(x+t)]  -  J  c-<j'2(t,cz)cdz 

=  2  [l2(x-t)  +  l°2(x+t)  +  BO(x-t)  -  EP(x+t)]  -  2n  J  f  a(t,y)dy 
Similarly, 

B(t,x)  =  |  [c^E^ctx-t))  -  c^EgWx+t))  +  c-'B°(c(x-i))  +  c1  B °(c(x-t))] 

cx  00 

2k  f  2k  f 

-  —  J  c-i  j  2(t,y)dy  +  —  J  c-i  j2(t,y)dy 

C  -00  0  cx 

=  ^  [  E°  (x-t)  -  E°2(x+t)  +  BP(x-t)  +  BP(x+t)] 

X  00 

-  2 n  |  C'1  j  2(t,cz)dz  +  2k  J  c-1  j  2(t,czx)dz 
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oo 


=.  |[E°(x-t)  -  E°2(x+t)  +  BP(x-t)  +  0>(x+t)] 


X 


-  2k 


I  j  2(t,y)dy  +  2k  . 


and  so  (RVM~)  with  c  =  1  is  satisfied. 


]  2(t,y)dy 
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